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THE RADIATION FIELD IN A FLUID IN MOTION 
By L. H. THOMAS 
[Received 20 August 1930] 


1. Introduction. The interaction between radiation and fluid in 
internal relative motion has been treated from various points of view 
and with various objects by Jeans,* Eddington,+ Rosseland,t Vogt,§ 
and Milne.|| These authors have usually neglected explicitly or tacitly 
small terms arising from the finite velocity of light. It is not difficult 
to extend the treatment to include such terms exactly, and it then 
appears that the terms neglected were sometimes of the same order 
of magnitude as those retained. This arises mainly since the dif- 
ferences in the radiation-intensity J, referred to a fixed coordinate- 
system and referred to a coordinate-system moving with the matter, 
due to aberration, Doppler effect, and the transformation of the 
energy-density as a component of the stress-energy tensor, are all 
of the same order of magnitude. 

In the sequel equation-systems invariant under Lorentz trans- 
formations will be obtained which take these effects exactly into 
account. 


2. Notation. The material stress-energy tensor and 
equations of motion. Of three notations, 
(a) the general relativity notation with constant fundamental 





tensor, 
—c?, 0, 0, 0 
a 0, 1, 0, 0 
Iv~=) 0,0,1,0[’ 
0, 9, 0, 1 


> > > 


(6) the usual vector notation with (a.b) for a scalar product and 
I for the unit dyadic, 
(c) ordinary Cartesian notation, 
that one will be used which is, at any stage, most convenient. Nota- 
* Monthly Notices, 85 (1925), 917; 86 (1926), 328; 86 (1926), 444. 
+ Monthly Notices, 79 (1918), 13; Internal Constitution of the Stars, p. 90. 
{t Astrophysical Journal, 63, 342. 


§ Ast. Nach. 5545 (1928), 232. 
Quart. J. of Math. (Oxford) 1 (1930), 1; Monthly Notices, 89 (1929), 518. 
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240 L. H. THOMAS 
tion (a) is most useful for expressing general equations of motion, 
but is conveniently replaced by (b) for dealing with radiation- 
intensity, while (c) is convenient for special cases and some proofs. 
The ordinary material velocity will be denoted by 
(u,v, w) = V, 
so that 
(1, u,v, w) (1,v) 


V4, V#, V3, V3) = —— Spee SUR <a 
sii ) {1—(u?+-v?+w?)/c?}* = (1—v?/c?)t 


defines the contravariant velocity four-vector V+. Likewise 
(26, 2,6, 26) _ (0b, ab oF 28) _ 4 up 
0x,’ 02,’ tra’ OX, at ’ da’ dy’ ez Bi a 
where ¢ is an invariant, defines the covariant gradient of ¢, be ; 
Cty 
For a perfect fluid, that is if we neglect self-diffusion, conduction 
of heat, and viscosity, the state is specified by the velocity v, the 
temperature T, and the invariant particle-density po, i.e. the particle- 
density referred to a coordinate-system in which v is zero at the 
point considered; the pressure p and invariant energy per unit- 
number of particles U being functions of py and T. The mass-density 
referred to such a coordinate-system is then 


Poo = Po( 1+ U/c*), 
Pooo = Poo TP, c 


the particle-density and flow-vector is pyV#“ and the material stress- 


and if 


energy tensor is pog9V#V"+-pgt”. 
The five hydrodynamical equations are, then, 
(1) the equation of continuity or conservation of particles, 





(pV) = 0 (1) 
Ty 
or, if 
D _ on? o l 2 ee ee 
Dt éx,  (1—v?/c?)! & a a )* 
ba ; ane = Cisne (1.1) 
po Dt OXy at (1—v*/c*)* (1—v?/c?)# 


(2) the four equations of conservation of energy and momentum, 


“ (PoooV4V"+-pge”) = FY, (2) 
02 y 

















RADIATION IN A FLUID IN MOTION 241 
where Fe = (T /c?, F) = (T/c?, X, Y, Z), 
and 7, X, Y, Z are the rates of increase of the energy and momentum 
of the matter per unit-volume. 
If py = pooo/(1—v?/c?), equations (2) take the more familiar form 
0 p é é é 
wa ome ne pest ask ; ~s = T 2 
ail Z)+ geet spi Sas [c 
é é é 0 
pg hret ax (px P)+ 5 awe — x 
é é é 7 _ 
at piv Ox pyvu+ ay (pye?+-p)+ r= } 


7 é é é 
agree Oe ate ayer 5; (Pi +P) = Z, (2.1) 


or, more shortly, 


o 
at (3) +V. py = Tie? 


a) 
ahiv + V «tev + pl} =F (2.2)* 


If equations (2) are multiplied by —g,,V° and summed, we obtain, 
since g,,V’V" = —c?, 





: PoooV#—V4 P = —9,, FV", 

eo # 

or, subtracting c?+ U-+-p/p, times (1), and writing 7’, for the value 
of T in a coordinate-system in which v is zero at the point con- 


(3) 


c2 


sidered, 
D yy_PPpo_ 


yg a 

For ordinary body-forces (such as gravitation, apart from terms 
due to general relativity) 7’) = 0, and equation (3) is the one that 
gives most simply the change of temperature at a point, 

The above has been written out to illustrate the notation; while 
the terms depending on ¢ are not usually required for matter, corre- 
sponding terms for radiation may be important and it is convenient 
to have the material equations in a form invariant for Lorentz 
transformations. 

* See Eddington, Mathematical Theory of Relativity, 116, 121, where, how- 
ever, py is used with a different meaning, 

R2 








242 L. H. THOMAS 

3. The effect of a Lorentz transformation on radiation 
intensity; Schwarzschild’s equation. The intensity of radiation 
I, of frequency v at x,y,z at time ¢ in direction with direction- 
— 1 = (l,m, n) = (cos 8, sin 6 cos ¢, sin @ sin ¢) 
is a function of x,y,z,t,v,0,¢ such that the radiant energy of frequency 
between v and v-+dv moving in a direction within solid angle dw 
about 1 through area dS in time dt is in the limit 

I,(1.dS) didwdv. 


On the electromagnetic theory this description assumes that we have 
averaged over phases in some way; and the form is consistent with 


a particle theory. 
Now consider the Lorentz transformation 
x’ = (w+qt)/(1—q?/c?)! 

a 
— ff 
=i 


, 
, 


, 


+(gx/c*)}/(1—¢?/c*)*, 
in which a particle at rest in bie is moving with velocity (q, 0, 0) 
in (wiy'z'\t’). 
Since ct—x cos @ —y sin 6 = en/v 
becomes 
a, cos 0+q/c nel —q*/c*)* _ en(1—q?/c?)# 
l+qceos@/e “~ 1+ qcos6/c v(1-+-q cos @/c) 
we have for a ray with direction given by (9, ¢), (4’, ¢’), 
¢=¢ 
cos 0’ = (cos 6+-q/c)/(1+-q cos 6/c) 
sin &’ = sin 6(1—gq?/c*)*/(1-+-q cos 6/c) 
v’ = v(1+¢qcos 6/c)/(1—gq?/c?)!, 
equations giving aberration and Doppler effect.* Notice that 
(1+-q cos 8/c)(1—q cos 0’ /c) = (1—q?/c?). 
Since dw = sin 6 déd¢ = dcos 0 d¢, 
dw’ = dw(1—gq?/c*)/(1+-q cos 8/c)?. 
Further, 
dt’ = dt(1—q?/c*), (dx’ = 0), since t = (t’—qza’ /e?)/(1—q?/c?)* 


* Cf. Jeans, loc. cit. 2. 
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The number of quanta of the type considered passing through area 
dS perpendicular to the x-direction and moving so that dz’ = 0 is 
proportional to 
(Ii /v') dw’dt’ cos 6’ dSdv’, 
and this must be the same as 
(I,/v) dwdt cos 6 dSdv-+-(I,/v) dw (q/c) dtdSdv, 


the number which would have passed through dS, if it had moved 
so that dz = 0, together with those contained in a volume of base 
dS and height q dt, since (I,/cv) dwdv must be the number contained 
in unit volume, and since qdtdS is the volume referred to 2,y,z,t 
swept out by dS. 

Hence I; = 1,(1+ cos 6/c)3/(1—gq?/c*)!, (4) 
and I’, dw'dv' = {(1+-q cos 6/c)?/(1—q?/c?)} I, dwdv. (4.1) 

These results extend at once to the Lorentz transformation in any 


direction 


r= {rater +1044 eh) arata:| /a—atjey 


t’ = {t+-(q.r)/c?}/(1—q?/c?)}, 
where r = (2, y, 2). 
The transformation-formula for direction-cosines of a light-track is 


, 2) n2 ei 2/2 1. | / i 4 
’ = {i —a"/e*)+1—(1—at/e2)]' Pata / {1+-(1.q)/c} (5) 
while the factor (1-++-q cos 6/c)/(1—g?/c?)? becomes 
{1+ (1.q)/c}/(1-.q?/c?)!. 

It is clear from the formula of transformation of J, that there is 
no inconsistency between the particle-description of radiation and 
restricted relativity. 

Schwarzschild’s equation for the transmission of radiation may 
be written 

l: 
~ +l. VL, = —pk, I, + pe (6) 
Generally, if it were necessary to take account explicitly of scattering 
and fluorescence, we should also have a term of the form 
[ [ pe(v,v’, 1.1’) D,(1’) dv’dw’ 


on the right-hand side, and this could be treated in the same way. 

















244 L. H. THOMAS 
Now 
se tl'.V) =[0—atle1+d.ay/9i(25+0.¥)), 
c ot c ot 
so that 
—(pk,)'I), +-(pe,)’ = [{1+-(1..a)/e}?/(1—q?/c?)](— pk, I, +-pe,) 
(pk,)’ = pk,(1—q?/c*)!/{1-+-(1.q)/c} 
(pe,)’ = pe,{1+ (1.q)/c}?/(1—q?/c?). 
If, therefore, the values of pk, and pe,, referred to a coordinate- 


and 


system in which v is zero at the point considered, are pok,. and po¢,o, 
equation (6) becomes 
] 7 2/,9\1 
~T,+(.V)L, = —pokyf {1—(v . 1)/c}/(1—v?/c?)}* 7+ 
F 
+ po€vo(1—v?/c?)/{1—(v.1)/c}*. (6.1) 


So long as we neglect very small terms depending on the departure 
of the velocity-distribution at a point of the matter from the Max- 
wellian, we may suppose k,, and ¢,, to have their equilibrium values. 

4. The radiation stress-energy tensor. The rates of absorption 
of energy and momentum by the matter from the radiation are 


oP ire 
| { pk, I,—pe,} dvdw, - | {pk,I,—pe,}1 dvdw, 
C 
and it is verified at once that 
ire, ri, ‘ 
(<— | | {pk,I,—pe,} dvdw, ~ | {pk,I,—pe,}l dudes) 
c* s c * * é 


transforms as a contravariant four-vector. 
Likewise 
“ire ire 
_ itten = | Ll dvdw 
a oe 


= | | ZA dide, : [ [ LAL drde | 
a c * * - 


transforms as a contravariant tensor, and by (6), 


os [| I, dvdw +V.5 | [ 2.1 dvdw " I {pe,—pk,I,} dvdw 
ce | | 2 


dt c® J , 


LA | | Tldvdw +V. all TM dvdw = all {pe,—pk, 1} dvdw. 
auc” J y c c 
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Thus this tensor is the radiation stress-energy tensor and either 

side of these last equations may be added to the left-hand sides of 
equations (2.2) to replace the radiation terms in them. 

Moreover, if one diagonal term of a symmetrical tensor is known 

in all coordinate-systems, the tensor is determined, so we need only 


actually calculate a [ | I, dvdw and can infer the other components 
c 


from it. 

To obtain J,, we solve (6.1) by successive approximations, ob- 
taining 

(1—v2/c?)? eo _ 
sii {1—(v.1)/c}3 kyo 
_ (l-vJe# 1 ge ) (1—v3/e*! ce 
{1—(v.1)/c} pok,» \c at }{1—(v .1)/c}3 kyo ; 

where except under conditions of rapid variation* the two terms 
written down give a very good approximation 


The largest terms in 5 [ | I, dvdw are therefore given by 


1 aie (1—v4/9¥ 66 aydey 


{1—(v.1)/c}3 kyo 
which is conveniently evaluated in terms of v, and 1, as 


Ed [pees €0 
C3. 
or, if 


(1—v2/c*) kg 
Ey dvy = - E,, since [ow .1,) dw, = 0, iG 1)? dwy = 
‘ 7 


v0 


dvgdw, 


e 


4x 


@ lve  ° @|31—v3e 3 
from which the whole tensor has the form 
[ 4 1 4 Kw 
| 3c? 1— vic 3c ° 3c? 1—v?/c? 
4 Ev 4 E,vv 
40 oh ee ‘0 sip 


| 3c 1—v?/ C2 3c? 1—v?/c? 


= 1 
Ba (53) ey "+ 30"). 


1 1+-v2/3c? 1 1 
of? 


* See Jeans, Monthly Notices, 86 (1926), 574. 











246 L. H. THOMAS 
The terms to be added to (2) can at once be written down, 


é Mu L[/v + Lv 
2% o(saV* ao" ), 


Oa 
Ox, 
! 


and those to be added to (3) are 


4 ool (E,V+)— 1 Vek, 
3 OL y 3 OX, 


c 


° 2 Baga VAV' +59") 
OX, ; 


—GugV 


DE, _ 4 Ey Dpy 
Dt 3p, Dt 


on subtracting 4H,/p) times (1), so that equation (3) becomes, if 
terms of higher order are neglected, 
De: aie 4E,\D a ; 
Por 0+ — B—(£+-—9) = 7 (3.2) 
Dt Dt Po 3 po Dt 
where 7', now contains no terms arising from either the ordinary 
temperature-radiation or from the gravitational field. 

To obtain the terms of next higher order, which include radiative 
transfer of heat and radiative viscosity, we proceed in the same way, 
making use of the formulae 

] 
— | (1l,.a)(1,.b) dw, = (a.b), 
47 | : 
LF 
— | (l,.a)(1,.b)(1,.c)(1y .d) dw, 
4. 
=} {(a.b)(c.d)+(a.c)(b.d)+(a.d)(b.c)} 


For higher terms further formulae are required, which are exactly 
analogous to the above; the next coefficient is 1/3.5.7 and all possible 
partitions of the six constant vectors into pairs occur. The corre- 
sponding integrals with an odd number of products are zero. 

The next term in J, is 


_ i— c2)jt 1 (= <4 
{i—(v. 1) |} pokye 


Supposing now that e,9/k,») is a function of vg and T only, so that 


a(e2) = sala) @ +58) 


ary 


v)) (l—v2/c?)# ¢ 
c at {l—(v. 1), c} kyo 
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c at 


we obtain 
(1—v2/e?)2 1 8 (e\ (12 
{l —(v. 1)/c}4 dive =e) (; at +(1. v) T+ 
ts ot A AM 1—(v.1)/ 
/c2) ( att. )(q ru). 


kyo) | Q—(v De} 


kyo 
£0 __ s¢{” 
kyo s #(:) ; 


30 a(S st) <0 5 
kg’ dv0\Kvo 


So _,, © (S0 
+1332 rez ) 


or, if we assume 
4 ( £0 
oT \kyo 


so that 


we have 
3 en la 2 (62) x 


 {l-(v 1) /c}# pokyo Ov 
) 


mt at +e. v))x— ( Saat (- at e. v))v} 3 e 
tases: Gat at ©: v))¥ }]: 


If we write 
1 @ ° 1f 2a (e 
0) dy, = - A (s) avo 

kg OT (= ) ” =| x vo "° 


so that « is the coefficient of opacity and is in general a function of 
‘ —y2/e2\t 
, i is 
(1+v.1,/c) 


and ppo,* the contribution to . | I, dvdw = — | I 
c c 


aa € c 
is, since [ 2a», = 2 E,, 
Jk, 4r 


v0 
1 0, {1-+( (v.1,)/c} 
- -~—+(l. 
C*xpy OT rod ie (1—v?/c?)? (es ait | v))r+ 
put (v .1,)/c} (1 ) ) 
T (1—v2/c2) |e’ (55 at (. de of 
oe fy (- z d 
If in this expression | is replaced by its value in terms of 1, and v 
t+-{1—(1—v?/e?)A (I, -Viiviiv+vie 
1+(v.1,)/c 
* Rosseland, Monthly Notices, 84, 


T 


1,(1—v?/c?) 
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and the integrations are carried out, the result is 
1 @E,f1 l+v%/c? or, 1! 
c (1—v2/c2)! at | ¢ ( 
_T 1+4}v?/c? ( 
c3 (1—v?/c2)i \ 


C*xpy OT 


1 GE, 2/c ; 
c*xpy OT | (1—v?/c*) (1—v?/c?)! 7 


x (-+ (w.¥))2— wa ‘ (v.V))2— 


at 1—v?2/c2)! \at © 


ny | 


—~T+—|— +V. : 
(1—v?/c?)t ot 5c \et (1—v?/c?)! (1—v?/c?)? 


2/3c O. . O& (: 1 Vv v )- 


I (: 1 Vv 2T a 1 
15c \at (1—v2/c2)! * * * (1—v2/c2)!) 5c a (1—v?/c2)! 
from which the whole tensor has the form 


1 0, 


Cxpo oT 


pe Es 
i 


‘ines ee ae ae ee oT 
2) Ky vy} pate ES Oe gt#¥c?] T | _¢2} i 
éx, | 3° ox, 3 Gat 


Xs he “+ 


Pee Po ok ee 1 oO = 
+ = c2) elt : +—T- (J LVVy T)+-— ger - a Te 
3 6x, 5S On, 15 Ox, 


L « 
= pee c'ngHt : 
15 


C2 OX 


T 


C VY + ong = ve. (8) 
a 5 ° 


the largest terms of which are 


? . oe - Dr _ 2 wie ity v) . ovr 
C’Kpy OT 3 Dt oa ss 3 


i 1 ~ | 
+—c?Vv+—c?Vv | 
15 15 


(8.1) 
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The terms to be added to (2) and (3) follow immediately. In 
particular the principal term in ‘ is 
—V..—_V,, (9) 
i 
giving 
D D Dp c 
pec © Ul = ae 77. E,=T,. J 
PoTt + 7 Mo -(2 13 3); Dt * "rh o= Ts ites 
5. Comparison with previous results. Equation (3.2) above 
agrees with the result of Jeans as corrected by Vogt*; the 


result given by Milne} in which 5 E, is replaced by © E, and 
o 


— ; a by — ; =! ~ requires, neglecting small relativity correc- 
tions, the addition to the ‘radiative conduction’ terms (9) a term 
(v.V)£,+ £,(V.v) or (V.v£,); i.e. Jeans’s F must be c(3xp9)-? V Ep, 
(a), while Milne’s F must be vE,—c(3xp,)-! VE, (b). The value of F 
given by Milne,t however, is $ vE,—c(3xpy))-! Vp, (c). In fact (a) is 
the rate of flow of invariant energy (energy measured in a coordinate- 
system moving with the matter) in a coordinate-system moving with 
the matter, (b) is the rate of flow of invariant energy in a fixed co- 
ordinate-system, while (c) is the rate of flow of coordinate energy 
(the energy component of the energy-momentum vector) in a fixed 
coordinate-system. 

The terms in the stress in (8.1) depending on the temperature- 
gradient agree with those given by Milne and differ from those given 
by Jeans because, as in Milne, but not as in Jeans, T here denotes 
the material temperature. 

As regards radiative viscosity, however, the state of affairs is 
different. The principal terms in the stress-tensor representing visco- 
sity are 

_ 1 ,0My[ gu, 20, om au, yu, 

lickpy OT Ox Oy o& Cy Cx Cz (Ox | 
duo gu = ew , @ 
cy On: Cy Ox Oz cy && 

ow ov , ow ou , Ww 


ste )46|OUktm)6OCRtete 


oz Ox cz (Oy 





Loc. cit.; see also Jeans, Astronomy and Cosmogony, p. ian 
Quart. J. of Math. (Oxford), 1 (1930), 8. 
Monthly Notices, 89 (1929), 523, equation (28). 
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The cross-terms have twice the coefficients found by Jeans and Milne, 
and the diagonal terms twice the coefficients found by Milne, Jeans 
differing further on account of his different definition of tT. The 
corresponding term in the energy-density is likewise doubled. The 
error in Jeans’s calculations occurs in the equation before his approxi- 
mate equation (7).* The term 


2 : (du sin @ + dv sin 0 cos ¢ + dwsin @sin d) 
- 

should be doubled, the extra part arising from the change of co- 
ordinate-energy from one.coordinate-system to the other. That in 
Milne’s is the same; his equation (5) is wrong and the term 2(i. v)/c 
in (6) should be doubled; his equation (14) is not invariant for Lorentz 
transformations and is not equivalent to (6.1) above. 

6. On the inclusion of generation of energy in the equations. 
If generation of radiant energy at the expense of the number of 
particles in the matter is to be included in the equations, a term 
representing the disappearance of particles must be included in (1); 
as well as corresponding terms in (2), say —p G@/c? and —pp)GV"/c? 
on the right-hand sides of (1) and (2). Terms representing the absorp- 
tion of the radiation produced would have to be included also. This 
radiation, however, would not be part of the material temperature- 
radiation considered above, but would be treated separately by 
making assumptions suited to the physical theory of generation on 
which it was based. For instance, it might be assumed that G repre- 
sented a property of matter independent of temperature and pressure 
and irreversibly producing high-frequency radiation absorbed by the 
matter at a definite rate pp»K. The above method could be applied 
to find the corresponding terms in the equations, f po¢,9 dvy being 
replaced by po)/47; and the first-order terms to be added to (2) 
would be 7 1Q4/ 


é pepe Leo! 
ea, ¢ K \3c* V3? ) 


and to (3), in addition to the corresponding terms, 


ior 4, SER. 
ay rarer she 32K) 4 


on account of the extra terms in (1). If we keep in the largest of 


* Monthly Notices, 86 (1926), 446. + Ibid. 89 (1928), 520. 
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these terms only and insert the corresponding conduction-term, 
(3) becomes 

D D 1G p ,4H,,41 G\D 
ol +54 ( Bote g)—(2+5 a+ 5a RK) Dt 
1 oe a. | 

“ers (sem “ar 3cKp, ¢ x) 

In general, however, scattering should probably be taken explicitly 
into account. 

7. Summary. The equations of transfer of radiation in a fluid in 
motion are obtained in a form including terms of all orders in the 
ratio of the velocity of motion to the velocity of light, but the calcula- 
tions are carried only to the first-order terms, comprising radiative 
transfer of energy and momentum, in the coefficient of opacity. The 
results differ somewhat from those obtained previously by others. 
In particular the radiative viscosity is shown to be double that found 
by Jeans and Milne. 


(3.3)* 


* This equation is, of course, essentially identical with either that of Jeans 
and Vogt or that of Milne, provided that the conduction-terms are written 
correctly. 








J CONVERGENCE AND SUMMABILITY CRITERIA 
FOR FOURIER SERIES+ 


By J. J. GERGEN 
[Received 22 June 1930] 
PART I 


1. Introduction. The first object in Part I of this paper is to estab- 
lish the generalization, stated in Theorem I below, of the Lebesguet 
criterion for the convergence of a Fourier series. The second object 
is to discuss the relation of this test to the six commonly recognized 
criteria and certain of their generalizations. It is found that in this 
form the Lebesgue test strictly includes the others under considera- 
tion. Incidentally, it will be proved that the continuity condition in 
one of the alternative forms of Lebesgue’s original theorem is 
redundant. 

In Part II we extend Theorem I in two directions, obtaining in 
Theorem II a somewhat more general convergence criterion and also 
a criterion for Cesaro summability§ of order p. Theorem I is con- 
tained in Theorem II; but it seems best, even at the expense of 
repetition, to give a separate proof. Numerous details, which tend 
to obscure the arguments and which do not occur in the proof of I, 
arise in that of II. 


+ The author is indebted to Professor Hardy, not only for his encourage- 
ment and his interest in this work, but also for suggesting, throughout the 
preparation for publication, several possible simplifications in the proofs and 
generalizations of the material. 

t Lebesgue, 13. 

eo 
§ The series > wu, is summable, to sum s, by Cesaro means of order p, where 
n=0 
p> —1, if 


n 


im [tn > ae 
= (P) > = nln] =e 
A’? 


n> 2 
m=0 


ol ; PF (n-+p)(n—1+-p)...(1+p) mee 
where Aj?’ = 1, and A‘?) = ——£—"__— <<< Pl tor n > 0. Summability 
n! i 
of order zero is of course equivalent to convergence, and summability of order 
p implies summability of order p,, to the same sum, if p, > p. 
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To simplify the writing, we shall suppose once and for all that the 
Fourier series 


f(t) ~ dao+ > a,, cos nt 


n=1 
in question is that of an even function f(t) which is integrable in the 
Lebesgue sense and periodic with period 27. Moreover, we shall con- 
fine our attention to the behaviour of the series at the origin, that 
is, to the behaviour of the series 


hay + > a,. (1.1) 
n=1 
We write do(t) = o(t) = f()—s, 
t 
1 
and (t) = = t—u)’6(u) du 
4,0 re] y-Ap(u) 
for r>0, where s is a number whose value will be assigned when 
the occasion arises. We denote by (C,), (C¥), and (C) the conditionst 


$,(x) =e (2”), 
d(x) = [isd] dt = o(), 


0 
(C) the series (1.1) is summable, to sum s, by some Cesdro means. 

In order that (C) hold, it is necessary and sufficient, as Hardy 
and Littlewoodt have shown, that (C,) hold for some positive 
integer r. 

The conditions sufficient for the convergence of the series (1.1), we 
recall, are then in 
(D) Dini’s test:§ $(t)/t is integrable on the interval (0, 7); 

(J) Jordan’s test:\| f(t) is of bounded variation on some interval (0, &), 
where 0 << E<7;3 
(V) de la Vallée-Poussin’s test:+} the mean value 


1 
F(t) == | f(u) du 
| 


is of bounded variation on (0, €); 


+ The variable x is always positive, and the symbols o, O always refer to 
the limiting process 7 — +0. 

t Hardy and Littlewood, 5. § Dini, 3. 

| Jordan, 11. tt de la Vallée-Poussin, 20. 
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(Y) Young’s test:+ (Cy) and 


(Y’) {\agef(o}| = O(@), 
0 


where the integral denotes the total variation of t f(t) in the interval (0, x); 
(Yp) Pollard’s generalization of (Y):t (C,) and (Y th 
(Yur) Hardy’s and Littlewood’s generalization of (Y C) and (Y’); 
(HL) Hardy’s and Littlewood’s test:|| (C) and 
¢ 
(HL’) | |A,f(t)|? dt = O(z) 
0 
for some p> 1, where A, f(t) = f(t+-x)—f(d) 
(L,) Lebesgue’s test: (Cf) and 
é 


\A..f 
(L4) [ Pat! a — o(a); 


ry 
x 


(L,) Lebesgue’s alternative to (L,): (Cf) and 


é 
(Lg) | a,{*)] dt = 0(1); 


(Lp) Pollard’s generalization of (L.):t* (Cy) and 
g 
(Lp) lim lim | a," 


ko 2->+06 
kx 


In (J) the sum of the series is f(+-0), in (V) it is F(+0), and in the 
remaining tesis s. 
The relations between these tests have been considered by several 


authors.{{ The following diagram indicates, with reservations enu- 


(t)| dt = 0. 
t | 


+ Young, 22. In 23 Young states that (Y’) and (C,), for any 7, are sufficient 
for convergence. No proof for the case = 1 seems, however, to have appeared 
until Pollard’s paper 14, and no proof for the general case until Hardy’s and 
Littlewood’s paper, 6. 

t Pollard, 14. 

§ Hardy and Littlewood, 6. See also Littauer, 25. Littauer proves (Yy,) 
by a method essentially different from that used by Hardy and Littlewood, 

|| Hardy and Littlewood, 7. The theorem was originally proved with € = 1. 
It is, however, a simple matter to deduce from this the test stated above. 

tt Pollard, 14. 

tt In the notation of the diagram, Pollard shows in 14 that (Y’) implies 











CONVERGENCE AND SUMMABILITY OF FOURIER SERIES 255 
merated below, all the implications which are known, so as far 
the author is aware. In this diagram, the notation L,=1L,+C* 
signifies the equivalence of the condition represented by the letter 
on the left to the two conditions represented by the letters on the 
right. A directed line running from one letter to another indicates 
that the condition represented by the former implies that represented 
by the latter. The full lines show the relations between the tests, and 














the broken lines the relations between the characteristic parts of the 
tests. Identical conditions are not connected, nor are the relations 


(Co) 3 (CT) 3 (Cy) 3 (C)T 
indicated. In the relations involving (J), s is understood to have the 


value f(+-0), and, in those involving (V), the value F(-+-0). 
The phrase ‘ characteristic parts of the tests’ used above demands 


(L}), Hardy and Littlewood show in 7 that (J) implies (HL), and in 5 that 
(Y’) implies (C) = (C,' 

so that (Yp) = (Yur) D (Lp). 

It is trivial that (Y)>(Y,), Gg) >s) (0,)> (CMe) 

The remaining relations, together with references to the previous literature, 
are given by Hardy in 4. He considers only the case in which f(t) is continuous 
at the origin ; but it is a simple matter to deduce with his methods the implica- 
tions shown. 

In a summary of this kind, it should perhaps be pointed out again that, for 
reasons often stated, the importance of any test is not less because the test is 
a particular case of, or equivalent to, another. In particular, the test (Yy,) 
deserves comment in this respect. The degree of generality of (C) over (C,) 
certainly places (Yy,) in a different class from (Ly), and gives it, because of 
the economy in continuity conditions, an importance which is quite inde- 
pendent of its relation to (Y>,). + Here ‘>’, as usual, implies ‘implies’. 

3695 s 
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a word of explanation. Most of the tests contain two conditions, 
one of which is continuity, or some generalization of continuity such 
as (Cf), (C,), or (C), and the other a condition more individual to the 
test in question. We refer to the two conditions generally as ‘the con- 
tinuity condition’ and ‘the characteristic condition’: thus in (Y) the 
continuity condition is (Cy) and the characteristic condition is (Y’). 
An examination of the diagram reveals that several questions per- 
taining both to the continuity conditions and to the characteristic 
conditions in the various forms of the Lebesgue test are unanswered. 
In particular, one may ask whether (C¥) may be replaced by (C;,) 
in (L,), or, more generally, (C,) by (C) and (Li) by 


(ip) lim lim ‘dt =0 


kw 2->+0 0 
ka 


= 
*|AS(O) 
t 


in (Ly).¢ The answers to these questions are found to be in the 
affirmative in the first part of 

THEOREM I. Jf (LR) and (C) hold, the series (1.1) converges, to 
sum 8. Furthermore, if (Lig) holds, then 

(C) = (Cj). (1.2) 

The second part of this theorem is of course a generalization of 
Hardy’s and Littlewood’s{ result that (1.2) is a consequence of either 
(Y’) or (HL’). 

To show that the conditions of the theorem, viz. 

(Lp) = (Lz) +(C), 

are implied by those of the tests enumerated, we need evidently only 
prove that (HL’) > (Lz) and (Li) > (Li). 
The first of these implications has, however, already been observed 
by Hardy and Littlewood;§ the second we shall prove in section 4.| 
This proof is effected by first showing that (Lp) implies 


x 


[ |(2)| dt = O(a); (1.3) 


ny 


0 

+ In 7 Hardy and Littlewood discuss these questions to some extent with 
reference to their test (HL). They do not arrive at any definite conclusion 
but express doubt that any of these changes can be made. 

t Hardy and Littlewood, 6 and 7. 

§ Hardy and Littlewood, 7. Their proof is given also in section 6 of this paper. 

|| It might be noted that (L},) does not imply (Lj). Witness the example: 
f(t) = logt. 
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and it is by a slight change in this argument that we deduce in 
section 5 that (Lj) > (C¥). 

Incorporating these results, and also the trivial relation 

(14) > (LR), 
into the previous diagram, we obtain the one following. It is interest- 
ing to observe that not only is (LR) a consequence of the conditions 


D J 





Y’+Cy 





= Y'+C, 








Y 
HL = HI’+C 
j ' 





of every other test, but so also is (LR) a consequence of every other 
characteristic condition. 

Perhaps a further word should be added in regard to the proof 
of I. First, the proof is of an essentially different character from that 
of Hardy’s and Littlewood’s proof of (HL), in that it involves neither 
the theory of functions of a complex variable nor negative orders of 
summability. Secondly, it is in part of the same character as the 
usual proof of (L,). The lemmas of sections 2.2, 2.3, and 2.4, in which 
it is shown that (LR) and (C) imply (C,), are of course without 
analogue in the proof of Lebesgue’s theorem, but the proof that 
(Li) and (C,) imply convergence (see sections 2.5, 2.6, and 3) follows 
the latter in principle, if not in detail, except at one or two points. 

2. Lemmas for Theorem I. 2.1. The proof of Theorem I rests 
on the following lemmas, some of which are known. In these lemmas 
and throughout the rest of the paper, AY%”®(t) represents the mth 
difference AmO(t) = Ss (—1)™49(™)@(t-++vx) 


v=0 


s2 


- 
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of the function ®, and 


: ae 
(2, k) = = ad — A, ~ 


ny 


ke ke 
We always suppose k > 1, and write 
a(x, k) = 6(1), 
for any function « of x and k, when, and only when, 
lim lim |a(a, k)| = 0. 
ko 2->+0 
2.2. Lemma 1. If (C) holds, then (C,) holds for some integer r. 
This is one of the pi i of Hardy and Littlewood cited in the 


Introduction. 
2.3. Lemma 2. If kx <v, then 
: ire : 
|AT+D¢, (t)| dt < (v-+-ra)ar [ —— dt (2.31) 
ka kz 


for every integer r> 0. 

This relation is plainly true when r is zero. Moreover, assuming 
it true when r is replaced by r—1, we have 
Pio t+vx 


fingng,co| dt = f Sayers) [ 4,-aw) de| ae 





kx ke * 0 
% » t+(v+1 
={/>(-1)"%) | $,-1(w) du| de 
kr |¥=9 t+ va 
v | tx (2 
=| | | 2 (1) )b,a(u-+vr) | du dt 
ke | t — 
v t+ax 
< | dt Jl AMd,_1(u)| du 
ket 
(k r1)x p 
= | (u—kx)|A%¢d,_,(u)| du+-ax [ia Wdh._s(u)| du + 
ke (k+ 1x 
U+zxr 
— J ( v—u+2)|A%d,_,(u)| du 
vw+2 U+rxr 
+ A blu)! 
<2 | |A®d,_,(u)| du < (v-+ra)a” | : — du. 
ka ka 


The lemma then follows by induction. 
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2.4. Lemma 3. If (LR) holds, and (C,) holds for some Ni r> 0, 
then $,(x) = 0(2). 
It is sufficient to prove that if r> 1 and 
r+1(%) = 0(2"**), (2.41) 
then ¢,(x) = 0(2”). 
We first observe that di de 
ha $,(a) = O(a, k)+ { dt J Av, _(u) du, (2.42) 
0 kt 
where h = (k+r)- and (2, k) is of the form 
D(x »k) = t dy sa( hkx) \tbd,safh(k -+1) )a}+.. +b,+1(2), 
the ’s being independent of x. In fact, 


hr x—rt ha ax—rt 


[at [ AP, s(w) du= fae f | 5 (—1*"(p)4,-a(-+¥0)} du 


0 kt 


he z+(v—r)t 


my fdt [  by-s(u) du 


0 (k+ v)t 


@)—-47)| dt +3 | [$,{0-+(v—r)}—d{(e-+v)0}] dt 


x hik+v)x 


he geoboae+S% f dat’ y f sae 


« hik+ pla 0 
= herd, (a)+ S bby .a{lk+v)2} = hag, (2) (a,b). 


Now, by Lemma 2, 
[han 2x—rt he 


| dt | Ay'd,_«(u uw <f @— tyr) dt fs 18db(u) 
10 kt 


0 


since kt does not exceed x—rt in the integral. ¥ . Hence, by (2.42), 
kt 


Sle) < < t, oe af ae b dt 


0 
ifa <€. But if (2.41) holds, the first function on the right is evidently 
o (1) for every k, while if (Lig) holds, the second is 6(1). Accordingly, 
IP, (a) |/a” = 0 (1) 
and thus, since ¢,(x)/x” is independent of i, ¢,(x) = 0(x"). This proves 
the lemma. 


har+1 
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2.5. Lemma 4. If (C,) holds, then 


a+kx 


i} BY) gin ™ dt = 0(1) 
t x 


for every a> 0 and k.F 

If a> 0, the lemma follows from the continuity of the integral 
f\P(é)| dt. If a =0, it becomes evident after an integration by parts. 
In fact, 


kz kx kx 
Pd(t) . at ae a { oy(t - 
$(¢) sin — dt = (kx) sin zk— — | Pr( leon dt + [ aa dt 
J ¢ x kx oF; % x J # x 
0 0 0 


dt = 0(1). 


é 
ee p(t+-2) = mt até 
B(x, &) = 2 lets as it = 6(1). 


We have 


2 
Vi 


ove Pt+2)’ 


a 
t(t-+-x)(t-+ 22) 


ig 


0v< 9i(%, t) rt at ee (t+ r) 


fort>0. Hence 


“ g 
B(x, k)| = $4(t-+2)g(e, sin sae... | $,(t-+2)g(2, t)eos dt + 
' Sin, 2: x 

kx 


rT | $(t-+-x)g, (2, t)sin ut dt 
"1 x 








ka 
g 
<— +e <5 
a J# k 


kx 
for kx < &, where M is an upper bound of {¢,(¢)|/¢ on (0, 2€). From 


this last inequality the lemma evidently follows. 


+t The proof of this lemma may be found in Pollard’s paper, 14. 
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3. Proof of Theorem I. If (C) and (LR) hold, then (C,) holds, 
by Lemmas 1 and 3. It suffices, then, for the justification of the 
complete theorem, to prove that (LR) and (C,) imply the convergence 
of the series, to sum s. 

Consider the integral 


a(x) = [Wi aa 


It is classical that the series converges to s if o(z)=o0(1). To show 
that this is the case, we set 


4 c+1)x k+2)x +x §&+ 
a(x, k) = {Ff a i + | | _2f - [ | nie dt, 
lo 0 0 é é 7 


and write é tia 


+22 
4o(x)—a(x, a j 4- | {fain dt 
kar (k+1)x (k+2)r ” 
f 
__ ot+z) 
t(t-+-a) —— zz) in™ dt 


kx 





FS j {a4 (ta) SAO sn! . 


t+-2x 


pas 2B(x, k)+y(x, k), Say. 
Now if (C,) holds, a(x, k) = 0(1), (3.1) 
for every k, by Lemma 4, and 
A(x, k) = 6(1), (3.2) 
by Lemma 5; while, if (Lg) holds, 
y(x, k) = 6(1) (3.3) 
because of the continuity of the integral {|4(¢)| dt and the fact that 
£+2x 
Iy(x, k)| < 2n(a, k)-+ 2 [ isiae 
g 
for kx < &. 
We conclude from (3.1), (3.2), and (3.3) that o(a) =o0(1); which 
completes the proof. 
4. Proof that (Lp) implies (LR). We first show that (Lp) im- 
plies (1.3). Denoting by p(x, &) the least upper bound of 


jt 
ky 
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for a fixed k, for 0<y<a, and, setting x,=2x{k/(k+1)}? for 





y=6,1...., we - 


ke ‘ 
[ is(t)| dt = > f i(t)| dt < 2 f _ dt 
0 peo key 44 a 
ws in eee ae ae TC) 
=k > 2,4 = 4 5 IPA)! dt 
2; (, J * kx, é ‘ 
=15] f Mla f Metsadla gf aoa 
hod 3 t+-2ya4 J 
yO Viet 43 katy +1 é ) 
i ftx al ld t 
<h> 2,4 a(x, k)+ | a dt} = 2(k+-1)- p(x, k)+- | ne dt- (4.1) 
v=o | : | : 


for kx <&é. But if (Lp) holds, p(a, &) = O(1) for some k. Hence 


(1.3) holds. 
Consider now 7(2, k). We have 


‘ é é 
| ASO yy < | A LPO)! at +o | HC+2)I oy 


J ¢ t } i(t+-2) 
kx ka k. r 
t+2 
($(t)) , 
A,{—} | dt l 
Pate H+ agers) | am, 
k 0 si 
2 | 
+ | cepa | Ml am 
ka 0 


r 


kx ‘ 
for kx < &, where M is an upper bound of : | d¢(u)| du on the interval 
x 
0 
(0, 2€). Hence, if (Lp) and (1.3) hold, so also does (Li). 

5. Proof that (L‘) implies (C7). 
This is an immediate consequence of (4.1). 

gta 


If (L5) holds, we have 


p(x, 2) = 0(1) and | - dt = 0(1). 
re é 
Hence J ie dt = xo(1); 
0 


and this is (C4). 
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PART II 


6. A criterion for summability. We obtain the generalization 
of Theorem I, mentioned in the introduction, by replacing (Lg) by 
the conditiont “ 

[AM F(t ‘ 
age, B) =a? [ 34 M at = 6 (1). (6.1) 
kx 

We prove 

THEOREM II.t Jf p> —1, and (6.1) and (C) hold, the series (1.1) is 
summable (C, p),§ to sum s. Furthermore, if —1 < p < 0 and (6.1) holds, 


then (C) =(C,), 
while if p> 0 and (6.1) holds, 
(C) = (E), 
where (E) represents the two canditions 
$,(x) = O(n), (6.2) 
$2(x) = 0(2*).|| (6.3) 


+ A generalization with mth differences of the Lebesgue criterion has been 
given by de la Vallée-Poussin. See 21, p. 150. 

{t The condition (6.1) is of course equivalent to that obtained by replacing 
the upper limit 7 by any number 0 < £ < 7, provided p> 0. This is no 
longer the case when p < 0; and it is, in fact, essential in the proof that the 
upper limit be 7. 

§ That is, summable by Cesaro means of order p. 

| When p> 0 the conditions (6.1) and (C) do not imply (C,). In fact, 
Hardy and Littlewood give in 10 an example of a function ¢(t) for which 


[ [pe] ae = 01m), lx) Aol@), —falr) = 0(@*); 
0 


and we prove below that = 


[Iso] ae = O(@) 
0 


implies (6.1) for every p > 0. However, it is worth noting in this connexion 


that (C) and w 
t 
a | A2fOl oy — o(1) 
tl+p 
x 
imply (C,) for every p, while (C) and (6.1), with p > —1, imply (Cy+ p+ g) for 
every § > 0. The first of these facts is a consequence of (8.26) ; and the second 
a consequence of the above theorem and Bosanquet’s result (see 1) that 


summability (C, p) implies (C,45+§)- 
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The proof depends largely upon the following lemma. 


Lemma 6. Suppose that —1<p<1, that S‘p(t) is the n-th Cesaro 
mean 


i/ 1 o 
S'P(t) = 4+ ag De es cos vt 
~n v=1 


@ 

of order p of the series 14. S cosvt, 
P a 
v=l1 


os ee 1 
that T'P(t) = 1 sin{(n- , +tp)i— spa} 
A'p) (2si $t)} +p 
and that R(t) = S'p(t)—T'P(p). 


Then, if 0<t < 3n/2, we have 
|S\P| < Mn, 


PS'P)| 
=| < Ma, 
dt? | 


where the M’s are independent of n and t. 

The first three of these inequalities are trivial. The last but 
one, (6.7), was first obtained in the general case —1<p<1 by 
Kogbetliantz,; by using Cauchy’s integral theorem. He shows that 

1 


__ sinpa (1+u)u"+e 
Rip'(t) = - ; we ——. da; 
2A ‘p) J (l—u)P(1—2u cos t-+-w*) 
0 


and from this equation both (6.7) and (6.8) follow. The proof of (6.7) 
given here is substantially the same as Szegé’st proof, which we 
reproduce in order to obtain (6.8). We follow this method chiefly 
to eliminate, so far as possible, the dependence of Theorem II on the 
theory of functions of a complex variable. The whole proof of (6.7) 
and (6.8) rests on the formula 1—e-“ = 2e(7-)sin 30, where 
i = v(—1), and on the elementary fact that, if p <1, the branch of 


+ Kogbetliantz, 12. For the formula in question see p. 277. 
t Szegé, 19. Szegé’s arguments are of the same character as those pre- 
viously used by Riesz, 16, in proving that 
MP)! < —pt—-\1+ p 
|S} <2 t f 
in the case 0 <p <1. 
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the function (1—z)!-* which is equal to +1 at the origin, is repre- 
sented everywhere in the circle |z| <1, with the exception of the 
point z = 1, by the convergent series 


1—2)-P = 14 5 AP-22, 6.9) 
> AY 


(v-+p—2)(v+p—3)... p(p—1) 


v! 





in which Ae-® = 


Two particular cases of Theorem II are worth noting: 
(HL,) Jf (HL’) holds, with =, and if (C) holds, then the series 
(1.1) is summable (C,—1/p+), to sum s, for every positive 5. More- 
over, if (HL’) holds, 
(C) = (C,). 


(HL,) If | \d(t)| dt = O(x) and (C) holds, then the series (1.1) is sum- 
0 


mable (C,8), to sum s, for every positive 5. Moreover, if | is(e)| dt = O(z), 
a 


(C) = (Cy). 


Both of these theorems are due to Hardy and Littlewood;} and 
although their prcofs as originally given are perhaps shorter and 
more illuminating than that of II, it is interesting that alternative 
proofs do exist. 

That (HL,) is a corollary of Theorem II is a consequence of the 


r) dt lq/a7 l/p 
i (jira 


ras A, f(t )| dt <a-} ral 
k 0 


inequalities 


pi/a+s8 
Ce 


lp 
M > 
<Al! fis svat 


0 


kx 


where p-!+q-!= 1 and M is independent of x and k. That (HL,) 


+ See Hardy and Littlewood, 7, for (HL,), and 8 and 10 for (HL,). They 
state (HL,) in 8 and prove it in 10. Pollard, 15, proves, by a method 
different from that used by Hardy and Littlewood, the particular case of 


xz 
(HL,), that (C,) and f |¢(t)| dt = O(x) imply summability (C, 1). 
0 
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is also a corollary of Theorem II is a consequence of Lemma 9 
below, and the fact that 


» {| 4, A (lt < Mae (4 SO! — yo Mae [ baad 
ti {ite } pte 
kx kz kx 
< Maue+ Mk-? = 6(1) 
if p> 0 and | |¢(t)| dt = O(a). 
7. Proof of Lemma 6. The proof of (6.4), (6.5), and (6.6) is 
immediate. We have 


since A?) — nP/T(1+p); (7.1) 
and similar inequalities hold for dS‘9)/dt, d?.8‘?/dt?. 

The proof of (6.7) and (6.8) is likewise immediate if p = 1, for 
gy — 1 {sin 3(n+-1)t)? 

“ 2(n+1)( sindt | 
sin{(n+1)t—}7} | 1 
A‘)(2sin $t)?  4(m+-1)(sin $t)?’ 

and so Ri) — {4(n+-1)(sin })?}-1, 

Consider, then, the values of p on the interval —1 <p <1. We let 


n 
on (z) = > AP”, 


n 
v=0 


and observe that S‘?) may be written in the form 
1 
S'P¢ 
nm ( )= ~ Alp) 
where i= ~+(—1) and R(U(2)} denotes the real part of U(z). We 
a riP)(z) = (1—z)!-P—(1—z)®olP(2), 
choosing that branch of the function (1—z)!+? which is equal to +1 
at the origin, and find, upon using the identity 
1—e-# — Q¢ei7-H2 gin Ut, 


~R[Seimfall (p) (e it) 4 giP) | (e -it)}), (7.2) 


that (7.2) may be written as 
A'?P)§\P) 


n nm 
eint = 1 ; , 
= R —— abe Rete + Wp (P)( pity y'P) _(e-it 
lame =r] ; 8(sin 4t)2 | Un ( ) n—1 )}] 
sin{(n+4+4p eden |g idee. 
F : gen (Ogle) §I- 
(Dain ty sip 8(sin 41)? [ rn (e—*) +r (e)}] 
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The first term on the right here is A‘?)7'?, by definition. Hence 
1 
8A‘?)(sin 3t)? 
Now, employing (6.9), we get 
re) : 
(2A‘ ~~). A‘P) + AlPa Prat z Ae a A(P))gnt2 > Ae-Sy 


v=nt+3 


R?(t) = pe , Re +My !P)(e—it) + ri?) (e-*)}). 


for all |z| << 1, except at the point z= 1. We thus deduce that for 
-9 
0<t<27 RP\t) = ViP(t) 
7 4A‘?)(sin $t)?’ 


nt+v+1 


ao 
where V‘?)(t) = A’P5Y+ > Am?) | cosvt. 
v=1 


This equation fs (6.7) immediately if we use (7.1). For we have 
|< { (m1 3 (ntytiyptl <=, 


) nt?" 


It also gives (6.8) if we differentiate after summing in Abel’s manner. 
In fact, 


Vip) — Aero 3 Ae 2) = >, AtsPasino-+0, 


in Mt 
and so 


wi FD beter 4 ES yn prp2y 


v=1 


n | a 


| dt | nr npe\n2-r2* nim 
< Mn-t 34+ Mn-4t-4. 
This completes the proof. 
8. Further lemmas for Theorem II. 8.1. In the rest of the 
paper y has the value y = n/(n-+43p+4), 


Hence 


jane _ MIVE)| xf l 1 


and N is a constant which may depend upon m and p and other 
parameters, but is independent of n and k. 


Lemna 7. If kx <v, then 


_ (m) | 
jaw | \idt< <a"(v-+rx)ir | — dt 


kx 


for every pair of integers r> 0 and m> 1. 
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The proof is in principle the same as that of Lemma 2. In the 
first place, the relation plainly holds when r is zero and m is arbitrary. 
In the second place, assuming it holds when r is replaced by r—1 
and m is arbitrary, we have 
v v t x 
| |Ar+md,(t)| dt = | | Ag+m-Pd,_1(w) du} dt 
kx ke 1 
t+2r 
a2 |Air+m-Dd,_.(u)| du 
kx 
v+r2 
Pr A@*D(t 
< xv'(v+ra)} p 2 (t) 
: fp 
kx 


dt. 


The lemma then follows by induction. 

8.2. Lemna 8. If (6.1) holds for some —1 <p < 0, and $,(x) = 0(x”) 
for some integer r> 0, then ¢,(x)=0(x). Further, if (6.1) holds for 
some p> 0 and ¢,(x) = O(2"), then ,(x) = O(x). 

We prove that (6.1), with —1<p<0, and 

$,1(%) = 0(x"**) 
imply d(x) = 0(2"), 
and that (6.1), with p> 0, and 
byaa(ee) = O(a") 
imply p(x) = O(2"). 
We first observe that 


he x—(m+r—1)t 
hh,xd,(x) = G(x, k)+- | dt | Aymtr-D¢g,_,(u) du, 

hye kt 
where h = (k+r+m—1)-, hy =(k+r-+m)-1, and G(x, k) is of the 
m—1 r+m—1 
Ss oa oe) oo 
> bdbristh(kt+v)a}+ d pdbypsrthy(k+v)2}, 

v=0 v=0 
the ’s being independent of x. This relation can, in fact, be verified 
in the same way that (2.42) was. 
We next observe that, because of (8.25) and Lemma 7, 


form ‘ 


Ge, k) 


hx 


,(x)| — |G(a,k)| , ae . ) 
ka. =: : t'-1-Py™(t, k) dt 8.26 
“UCT hh,at*} hh, . sided ( ») 


hyx 


for 0<2<7z. 
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Suppose, then, that (6.1), with —1 <p <0, and (8.21) hold. Under 
these circumstances, the first term on the right of (8.26) is o(1) for 
every k, while the second does not exceed the least upper bound of 
3"(a’, k) for 2’ on the interval 0 <2’ <hzx, and is therefore 0(1). 
Hence (8.22) holds. 
Suppose, on the other hand, that (6.1), with p > 0, and (8.23) hold. 
Then for some k, we have 7'7"(a, k) = O(1), and hence also 
hz 
sa Neoa )+O(z- | t'-1-p dt) = O(1). 
hyx 
This is (8.24), and — the proof. 
8.3. Lemma 9. If ¢,(2) th ) and ¢, (x) = 0(x") for some pair of 
numbers r,>r> 0, then es = o(a*+>) e every positive 5. 
This is a particular case of a ‘jens of M. Riesz.t 
8.4. Lemma 10. If r> 0 and ¢,(x) = 0(2”), the series (1.1) is sum- 
mable (C, r+-8), to sum s, for every, positive 5. 
This theorem for the case r= 0 is due to Riesz,{ for the case 
0 <r <1 to Hardy and Littlewood,§ for the case r= 1 to Young,]|| 
and for the case r > 1 to Bosanquet.t} All we need here is Bosanquet’s 
result. 


8.5. Lemma 1l. If 1<p<0 and (C,) holds, or if 0<p<1 and 
(E) holds, then 


(e+) 
lim lim | [ “$(t) Si?)(t) a0 (8.51) 
ape: @ 

for every v. | 

It is enough to prove the lemma with v= 0, since k+-v becomes 
infinite with k. We have, upon integrating by parts, 
ky | 
[ 40 S21 ae < da (ky)S (by) |+ fs b(t) Oe ae 
0 | 


= J(n,k)+K(n,k), 


t M. Riesz, 16. See p. 7. 

{ M. Riesz, 16 and 18. Chapman, 2, gives another proof of the same fact. 

§ Hardy and Littlewood, 9. 

|| Young, 24. Young’s result is a generalization of Lebesgue’s, 13, that 
(C,) implies summability (C, 2). 

++ Bosanquet, 1. Bosanquet’s result was stated but not proved by Hardy 
and Littlewood, 9. 
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Now, if —1<p<0 and (C,) holds, 

lim J(n, k) < N lim |nd,(ky)| = 0 

nn no 
= ky 

lim K(n,k) << N lim (n° ,(t)| dt) = 0, 

nD n->D 0 J 

by (6.4) and (6.5). Hence (8.51) is established in this case. 

On the other hand, if 0 <p <1 and (E) holds, then 


JI(n, k) < Nky{| TP (ky) | + | RiP (ky)|} 


< Nky{n-?(ky)-1-?-+-n-"(ky)-?} < Nk-p 

by (6.7), while 
: ds\P)}é 
lim K(n, k) <lim stn s | 


no nn dt a 0 


<N lim ‘ n2 \bo(ky)|-+-n | \p(t)| dt 
no 0 } 
From sie and (8.53) we deduce (8.51). 


8.6. Lemma 12. If 4,(x) = O(a), then 


| muy 


lim lim | i's f(t) R'?(t) dt| = 


by (6.5) and (6.6). 


kon n> | 
Kk+v)y 
for every p and v. 


We may suppose again that y= 0. We then have 


awtmy 
a 


awhpy 


” fla d Ri?) 
| d(t) RY (oat — t) RPE fie | b(t) at 


ky 


va 
for y small enough. The lemma then clearly follows. 
. Lemna 13. If $,(x) = O(a), then 


a—tmy 


: ; 1 . {tar pt “ 
ion iim | Ap d¢(t-+-vy)w sin ~ ) dt| = 0 (8.71) 
ky 
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for every m>v> 0, where 
m m—v v 
{sin }(t+vy)}'+" (sin }t)'*? = {sin }(t-+-my)}1+"" 
We first prove that 
My? dw| _ My? 
{3+p ” |dt| ~ tt+e 


w = w(y,t, p,v,m) =—_ 


(8.72), (8.73) 


for my <a and 0 <t<a—4my, where the /’s are independent of 
t and y. Consider the function 


nha m m—v v 
mw {sin £(t-+-va)}+e (sin 44)!+P ~ {sin }(t+-max)}+e 


as a function of x, for a fixed ¢ on the interval 0<t<7. Evidently 
w, vanishes with its first derivative 

F cos 4(t-+-vax) cos 4(t-+- max) 

x)= —}m(1 eee 

w1(#) 4mv{1-p) ta M(t+vx)}+" {sin }(t4+-mzx)}5+? 


at x =0. Moreover, 


v v(2+-p){cos $(t+-va)}? _ 


"(x)| = anv(1-+p 
<a dmv ) fein I (t+-va)}+P {sin }(t-+-va)3+P 


- m __ m(2+-p){cos }(t+-ma)}? 
{sin $(t+-ma)}'+P {sin }(t+-ma)}*+P 





< 2m*[{cosec (t-+-vax)}5+P+-{cosec }(t-+-mx)}*+?] < 4m(32/t)3+? 


for 0 <t<a—4mz and mx<7z. Accordingly, 


w(y,t, p,v,m)| = |w,(y)| 


This is (8.72). Similarly, the function 








mecos}(t-+-ve) _(m—v)cos at _—v cos $(¢-+-mz) 
{sin }(t-+vx)}2*? (sin }t)?*P ~— {sin 4(¢-+-mzx)}?+ |’ 


wz) = — a+e)) 

which for a fixed 2 = y is the derivative of w with regard to ¢, vanishes 

with its derivative at x = 0. Moreover, it is easily verified that 
w9(x)| < 3m5(3z/t)*+P 


for 0<t<a—4mz and mx <7. This plainly gives (8.73). 
3695 T 
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The proof of the lemma is now immediate. Integrating by parts, 
we obtain ‘ 


a— my 


1 j ‘tar ht 
ri | (t+ vy)w sin {= — sk dt 


ky 


] ta pt 7w—kmy 
try) sin (- ; )| = 


p) 
A n ky 


a—smy 
A 


$,(t- ey) Gsin(@ 4 —- dt — 
: dt 


7 dmy : 
7 ta | 
—- | d,(t-+-vy)w cos (= —FE dt 
y . y 2 
ky 
N 7 
. N " * dt : N 
< Ny?+P + Ny}+P < Ny?+P- 
oT t SY jap SY I 
ky 
for y small enough. This last function is annihilated by the operator 
lim lim; and the lemma is established. 
ko nx 
8.8. Lemma 14. If m’ ts an integer >m, and 7'%""(x, k) = 0(1), then 
gh (x ,k)= 7 o(1 i 
"he is sufficient to prove the lemma with m’=m-+1. We have, 
since f(t) is even and periodic, 


” |Agmf(t+-2) 


no t Na, k) < aP i - 


ny 


kx 


< My” (x, k)+a° | 7 


ny 


v7 


Arf (t) 


(227—t—x)4 


iat < Mnj(x, k) 


< M n(x, k)+-aP [ 


Tz 


for x small enough, where the M/’s are independent of 2 and k. The 
lemma evidently follows. 


9. Proof of Theorem II. If (C) holds, and if (6.1) is satisfied 
with p some number in the interval —1 <p< 0, then (C,) holds, by 
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Lemmas | and 8. Moreover, if (C) holds, and (6.1) is satisfied with 
some p> 0, then (E) holds, by Lemmas 1, 8, and 9. Accordingly, 
the proof reduces to showing that (C,) and (6.1), with —1<p<0, 
or (E) and (6.1), with p > 0, imply summability (C, p), to sum s. 

The case p > 1 offers no difficulty. If (E) holds, then 

Pii-+p) ~~ o (ait +?)), (9.1) 
by Lemma 9; and (9.1) implies summability (C, p), by Lemma 10. 
Suppose, then, that —1<p<1. The nth Cesaro mean s‘?) of order 
p of the series (1.1) is given by the formula 
(p) ae 2 St) d 
Sn eae p(t) n (t) t. 
0 


This we may write 


2m (ktv)y 2m aty—m)y 
= rer) [ sospma+ Ser | sORPH a+ 


v= 0 (ktvly 


aty—my 


Sem) fF sospqa+ ern f gorgat 


aty—my (k+v)y 


= a(n, k)+-B(n, k)+-y(n, k)+8(n, k), 
say. 
Now (C,) or (E) implies both 


lim lim |a(n, k)| = 0, 


k>~ n> 


and lim lim |A(n,k)| = 0, 
ko n> wo 
by Lemmas 11 and 12. 
Furthermore, since ¢(¢)S\?)(t) is even and periodic, we have 


y(n, ke) 
f smsprnat+ Sem f  semspwa 
v=m+ 


atly—my aty—my 


art(y—my 


MOSPQa— > er f BOSE a 


v=m+1 


awt(m—viy 


(9.4) 


T2 
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Hence it remains to consider 5(n,k). We write 


2m amy 
3(n,k) = 2) | d(t-toy)TP(t-+ry) dt 


ky 


7™—my 
(2m—1), te 
[ Ay2m-Dp(t-+y) sin( en Pom 


{sin 3(t+2my)}}+P 


aa ness 


ky 


m—mMy 


y 
(2m—1), / \ 
od [| = #0 sin (= Fj) dt +- 


(sin }t)!+? © a 


y 2 
ky 
2m—1 amy 
—ly, . A 
+ = ( ey. [ H{t-+ry)o sin (= —&) dt}, 
=m : y 2) 
y=] ki 
setting 


w= w(y,t, p,v, 2m) 
a 2m 2m — v v 
{sin K(t-+-vy)}? (sin it) e fsin (t+ 2my)}" +p” 
It is then seen from Lemma 13 that: 


lim lim |8(n,k)| << N lim lim 7"-"y, k), 


k->~ n> k>o nx 


if either (C,) or (E) holds; and further, from Lemma 14, that 


lim lim |8(n,k)| = 0 (9.5) 
kw n> 
if (6.1) likewise holds. 
The theorem now follows from (9.2), (9.3), (9.4), and (9.5). 
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A NEW WARING’S PROBLEM WITH 
SQUARES OF LINEAR FORMS 
By L. J. MORDELL 
[Received 2 September 1930] 

1. Ler f(X, Y) =aX*42hXY+bY? (1) 
where a > 0, h, b are given integers and 

A=ab—h?>0, (1.1) 
so that the quadratic form f(X, Y) is never negative for real values 
of X, Y. The new problem, which apparently has neither been men- 
tioned nor considered by other writers, is to find the least value of 
n independent of a, h, b, say N, such that integers 


a,, b, (r = 1, 2,..., 2) 


exist for which f(X, Y) = > (a, X+0,¥)? (1.2) 
1 


identically in X, Y. ‘ 

For one variable when } = h = 0 and also when A= 0, N= 4 as 
the result is simply Bachet’s theorem, first proved by Lagrange, that 
every positive integer is the sum of the squares of four integers. 

It is obvious that N > 5, on considering the case 

X?24-7Y2=— X24 Y24 Y24 Y244Y?, 

as then a, = 1, a,=a,=...=0, 6,=0 and 

7 = 63+63+62+6?-+... 
is not possible with less than four squares. In a paper to appear 
in the Mathematische Zeitschrift,* I proved that N<6 if h=0, 
N <8ifh+~0. The proof depended on the classic elementary theory 
of numbers, and was obviously not a method for finding the value 
of N. The difficulty was to find any non-trivial upper bound for NV, 
the value N < 12 being obvious on taking f(X, Y) to be a reduced 
form, i.e. a >b>2)\h|. For 

f(X, Y) = (a—Jh|)X2+ |h|(X+Y)?+ 6—|h|)¥?, 

and four squares suffice for each of these three terms. 

I also investigated there a different problem, namely, for given 
n = 2,3,4,5...in (1.2), to find the necessary and sufficient con- 

* ‘On the representation of a binary quadratic form as a sum of squares 
of linear forms.’ 











A NEW WARING’S PROBLEM 277 
ditions for the existence of rational values of a,, b,. It was shown 
that when n = 2, integral values of a,, b, existed if rational ones did; 
and similarly when »=3 provided [a,h,b]= 1, i.e. the greatest 
common divisor of a, h, b is one. The condition for rational solutions 
when n = 4 was shown to be that 

A # 4°(80-+7) (1.3) 
where p > 0, o > O are integers, i.e. that A can be expressed as a sum of 
three integer squares. It had already been shown by Landau* that 
rational solutions always exist when n = 5 and so of course if n > 5. 

This suggests that integer values of a,, b, exist when n = 4, if and 
only if (1.3) holds, and that they are always possible if n = 5. I have 
now found a proof depending curiously enough upon a type of 
theorem first proved by met some thirteen or fourteen years ago and 
dealing with the derivation of the representations of a number as 
a sum of four squares from those of the factors of the number. 

The proof of the subsidiary theorem is given in § 2; the case n = 4 
in § 3, and n= 5 in § 4. Section 5 is an application of § 3 concerned 
with the conditions under which numbers can be represented by 
certain ternary quadratic forms, while § 6 is a generalization of § 4. 

2. It is well known that if / is a positive integer, Jacobi has 
proved that the number of solutions of 

l=ajt+yit+a+ (2) 
in integers 21, Yj, 2, 4, 1S AD d, 


where the sum refers to all the odd divisors d, of 1, and A, = 8 or 24 
according as / is odd or even. 


If m = x3+y3+23+43, (2.1) 
where , 2», Yo, 22, tg are integers, then integer solutions of 
lm = x?+-y?+-2?+-#? (2.2) 


are given by Euler’s formula 
© = XyXg—Y1 Yo— 2%. —hfe 
Y = XyYoT Yet Aly—h22 (2.3) 
% = UyZq—YyloT2y%o +b Yo [ 
t= XyloTYi%2— 21 Yo Fh%. 


* ‘Landau, ‘Uber die Zerlegung definiter Funktionen in Quadrate’: Archiv 
der Math. und Physik, (3), 7 (1904), 21. 

+ ‘On the solutions of 2?+y?+2?+# = 4mm,’: Messenger of Math., 67 
(1918), 142-4. 
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The subsidiary theorem is that if [/, m]= 1 or - 2, all the solutions 
of (2.2) are given by (2.3), each solution occurring 8 or 24 times 
according as [1, m] = 1 or 2. The method of proof is the same as for 
my old result. Thus a unique correspondence can be established 
between a solution 2, y,, 2,, #4, of (2) and the linear substitution of 
determinant / in two variables. 


g =| XY, 24+, 
Sb, = . ° ? 
27h, %—WY; 


where the substitution i ; 


replaces X, Y rX+BY, yX-+6Y respectively. 


aes el . y Lot v4 - it. } 
Similarly for S,= ( at ey Zar ), 
—Z,+it,, Xy—Ye 


a xt+iy, z+at | 
whence SS, = ( ; Yy < ~ S. 
Z —2z+u, r—vy 
since from (2.3) 
oa ry (& + 1Y4)(%_+ Ya) + (2% +0ty)(—Z+ tte) 
+ it = (%+ ty) (Zit) + (% + 0ty)(%o— typ). 

Suppose now 2}, yj,... and x3, y},... are other solutions of (2), 
(2.1) giving new substitutions S;, S; and new values of gz, y,... in 
(2.3), say x’, y’,..... We wish to investigate when 

, , , 

“=, y =yY, Zz =2, t 

i.e. SS. at Sy 9 or S) 18 Si = SiSz! 
where S;" is the substitution inverse to S,, etc. Now w S'S, is a sub- 
stitution of determinant unity whose coefficients are rational numbers 
with a possible denominator /. So the denominators in S3Sz! are 
possibly m. Since [/, m] = 1 or 2, the denominator can only be 1 or 2. 

Suppose first that [/, m]—1. Then 


: E+”, 
1s, = | ef _ ] 
6-7 er, 
where £, », ¢, 7 are integers and 
tnt Etta. 
Hence = +1, A -0, (=0, r=0, ete., giving eight solutions, 
i.e. a value of x, y, z, t occurring once in (2.3) occurs exactly eight 
times. But there are A,> d,, ,,> d,, sets of values of x,, y;,... and 


Lg, Yo,-.., and so there are 
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BAA > > Im =Am > dim 
different sets of values of x, y,z,t. For, since 1 and m have no odd 
factors in common, 
X% > dn = dim: 

Also either Lat — ‘al 
or A=8, A,y=—24, Am=2%4, 
or A= 24, An=8, Nm = 24, 
and so AAm = 8A: 

Suppose next [/, m|= 2. Then 

S\18, = (_ ne in, 3 
30+ 37, 
where é, n, ¢, 7 are integers and 
&+7?+ +7? = 4. 

Hence wet gut, fin Hi, 

giving eight solutions, and 
@=1, y=1 @=1, =I, 
giving sixteen more solutions, e.g. 
é=1, y=], f=], r=1. 
These really give integer solutions for x}, y},..., for then 
wb tyy = 3(%,+0y,)(1+-2)+ }(2, +0, )(1+2), 
etc., where from (2), since / is even, 
2,+y,+2,+t, = 0 (mod 2). 

Hence a value of x, y, z, t oecurring once in (2.3) occurs twenty- 

four times. The same result as before follows since here 
A, =A, =A,,, = 84, J AAn = Am: 
3. It will now be proved that he idies of a,, b,, (r = 1, 2,3, 4) 


exist such that ‘ 
aX?+2hXY+bY?= > (a,X+b,Y)? (3) 
r=1 


identically in X, Y, ie 


4 4 
a= > a, b= > 8, at (3.1) 
r=1 


r=1 r=1 
if and only ifa>0, O<A=ab—h* # 4°(80+-7), (3.2) 
where p > 0, o > 0 are integers. 

The condition (3.2) is necessary. For we have identically, as a 


variant of (2.3), 
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aaa ae 
> a7 > bF = d ¢F, 
r=1 r=1 r=1 

where ¢, = a,b,+a,b,+a,b,+ a,b, =h, 


Cy = 4,b,—ayb, —azb,+a4bg, etc. 
and so A must be a sum of three squares. 

The proof of the sufficiency of (3.2) is simplified by some pre- 
liminary considerations. If (3) is proved for any quadratic form (1), 
it clearly holds for any form derived from (1) by a linear substitution 
with integer coefficients and any determinant. We may also suppose 
that the form (1) is properly primitive, ie. [a, h, b]=1. For if 
[a, h, b} =d > 0, then d is a sum of squares of four integers. Hence, 
if (3) holds for the form 


a h b 
ake _ xX Y4-_ 72 
a i ae. 


Euler’s identity shows that (3) holds for the form (1). This does not 
affect (3.2). 

Next we may suppose [a, b]= 1 or 2 by effecting, if need be, a 
linear transformation of determinant unity. 

The condition (3.2) states that A can be expressed as a sum of 
three squares or, say, ; s 1, = 
“" ab=h?+ U2+V2+ W?. 

Hence, from § 2, there must exist integer values of x,, y,, 2, ¢, and 
Le, Yo; 2%, t, for which 


a=ax+yi+22+#, b= 23+y3+23+4, 
while h is one of x, y, z, tin (2.3). Take, for example 
h = © YotYi%o+2to—h2e, 
then 
I(X, VY) = (@X+yo¥ P+ (YX +a2V P+ (qX+t¥)P+(hX—z¥). 
Hence the proof of the theorem, which is clearly independent of 
which particular value in (2.3) we take for h. 

It is well to emphasize that it has been proved that (3.1) is satisfied 
by integer values of a,, b, if and only if A 4 4°(80+-7), ie. the con- 
dition [a, b] = 1 or 2 used in the proof, is not required. 

4. It is now not difficult to prove the theorem that every definite 
quadratic form (1) can be expressed as a sum of five squares of linear 
forms with integer coefficients. It suffices to assume that 

A = 4(80-+7), (4) 
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as otherwise four squares would suffice. The theorem is true for 
A=7, 15. For 
X?47Y? = X24 Y24 Y24 y244y2, 
2X?4+2XY+4Y?2= X24+(X+4+ Y)2+ Y24 y24 y2, 
X24 15Y?2= X74 Y24+ Y244V2+9Y2, 
3X24 5Y?= X24 X24 X24 Y244Y?, 
2X2+2XY+8Y?2=(X+2Y)?4+(X—Y)?+ Y24 y24+ y2, 
4X21 2XVY+4Y?= X24(X+ Y)?+(X+ Y)?+(X—Y)24+ Y?2. 
Also any form with A equal to 7 or 15 is equivalent to one of the six 
left-hand forms. We may assume, then, that A > 23. On effecting 
a linear substitution of determinant unity and writing —y for y if 


need be, we may suppose that the form (1) is reduced and that 


tit hi b>a>2h, a<2WA7B. (4.1) 


Consider first the case p= 0 so that A= 8c0+7. Take one of the 
five squares to be q?Y*, where g is an unknown integer such that 
la) $ 
ei <>. So aX?24 2hX Y+(b—q?)¥? 
must be a sum of four squares. Hence from § 3 

0 <A, = A—aq? 4 4\(8u+7), 

say, where A > 0, » > 0 are integers. If 

a= 1, 2, 4, 5, 6 (mod 8), 
take q = 1 and then 

A, = 6, 5, 3, 2, 1 (mod 8). 
This requires A—a > 0 which is satisfied as a < 2vA/3. If 

a = 3, 7 (mod 8), 

take q = 2 and then A, = 3, 3 (mod 8). 


This requires A > 4a, and it suffices if 8VA/3 < A which is true as 
A> 23. 
We must still discuss a = 0 (mod 8) so that a>8. 
From ab = A+A?, 
b<jA+M/a, 
< }A+a/4, 
< 4A+4VA/3 < }A, 
since A>6. Hence, unless b= 0 (mod 8), we can take p*z? for one 
of the five squares where p = 1 or 2. 
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We must next discuss the case 
A=7 (mod 8), a= 0 (mod 8), b = 0 (mod 8), 
so that h is odd. We take one of the five squares to be r?(X+ Y)?. 


Hence (a—r?)X?+ 2(h—r?)X Y +-(b—1r?) Y? 
is the sum of four squares, and so r? <a and 
0< A, = (a—r?)(b—r?)—(h—??)?. 

= A—r*(a+b—2h) A 44(8u+-7). 
Now r? <a, 0< A, if r?<h, for then 

A, > ab—h*—h(a+b—2h) 

> (a—h)(b—h) > 0 

from (4.1). 


Also A, = 7+ 2hr? (mod 8) 


= 7+ 27? (mod 8). 


We take r=1. This requires 1 <h, which is true since h is odd. 
Hence the theorem is proved for A = 8c0+-7. 

Suppose next A = 4°(80+-7) (p > 0). By a theorem of Lipschitz,* 
every properly primitive quadratic form of determinant Dp*, where 
p is a prime, (and hence of determinant Dp?) can be derived from 
a properly primitive form of determinant D by a substitution with 
integer coefficients and determinant p (or p* in the second case). 
Hence it suffices to prove our theorem for the improperly primitive 
forms of determinant A, i.e. those with [a, 2h, b] = 2. But then A is 


even since A = ab—h?, and we can write 


aX?+2hX Y+bY? = 2[}aX*+ 2(4h)X ¥+ OY). 


The determinant of the form in brackets is }A. Hence, step by step, 
we are brought to the case p= 0. Hence the theorem is proved and 
N= 5in $1. 

A similar problem arises for the general positive definite quadratic 
form in s variables. In the Zeitschrift paper, it has been proved that 
the form can be expressed as a sum of squares of s-++-3 linear forms 
with rational coefficients. The coefficients are integers when s = 1 
or 2, and it would be interesting to know if this is true in general. 

5. We consider now the application of §3 to the representation 
of numbers by definite ternary quadratic forms. In the following 


* Matthews, Theory of Numbers, pp. 159-62. 
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scheme, the right-hand expressions, wherein p > 0, o > 0 are integers, 
denote those numbers and only those numbers which cannot be 
represented by the corresponding left-hand quadratic forms for 
integer values of x, y, z. 


xr y?+ 22, 4?(80-+-7) (5) 


x+y? + 22%, 
x+y?+ 32?, 
a? + 2y?+- 22%, 
x? 2y?-+-32?, 
22+ 2y?+ 422, 
x? 2y?+- 522, 


220+1(8g-1-7) 
3*11(3p +2) 
4°(80-+-7) 

220+1(8-4 5) 
220 +1(8g-1-7) 


5%+1(5-42) 


(5.1) 
(5.2) 
(5.3) 
(5.4) 
(5.5) 
(5.6) 


x2 y?+ 422, 4°(80+-7), 86-3 (5.7) 
x?-+ y?+ 52”, 4°(8a0-+-3) (5.8) 
x?+ y?+ 62, 37P+1(30-+-1) (5.9) 
a2 y2+ 822, 2%+1(8g-1-7), 1604-6, 40-++3, (5.10) 
x? 2y?+- 622, 4°(80-+-5) (5.11) 
a24 Dy2+1 822, 4°(80-+7), 86-+5 (5.12) 
x?+ 3y?+ 322, 9°(30-++-2) (5.13) 
x? + 5y?+ 52, 4°(80-+-7), 50-+2 (5.14) 
2a?+- 3y?+ 622, 4°(80+-7), 30+-1 (5.15) 
Qx2+ By24322,  — 9°(30-+-1). (5.16) 

Results of this kind are particular cases in the general arithmetical 
theory of ternary quadratic forms due to Gauss, Eisenstein, and 
Smith. This gives necessary and sufficient conditions for the repre- 
sentation of a given number by a genus of forms, i.e. a certain finite 
number of classes of forms. If there is only one class in the genus, 
the method gives at once the required conditions for the representa- 
tion of numbers by the class. 

Some of the results tabulated are due to Gauss, Cauchy, Lebesgue, 
Dirichlet, and Eisenstein, who very likely knew most of them. Exclud- 
ing (5.13) noted by Dickson and (5.14), (5.15), (5.16) noted by me, 
these results have been gathered together by Ramanujan,* who most 
likely found some of them empirically. They are not all independent. 
Thus (5.3), (5.7), (5.12) are obvious deductions from (5), while (5.5), 


* Ramanujan, ‘On certain quaternary forms’: Collected Mathematical Works, 
1927, pp. 171, 172. See also my remarks in the Appendix pp. 341, 342. 
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0) follow easily from (5.1). A proof of (5.2), (5.4), (5.6), (5.9), 
3) has been given by Dickson* and is based on the general theory. 
he theorem of §3 states that integer solutionsa,, b,, (r = 1, 2,3, 4) of 
4 1 4 
a= Se, k= 5 a)., b= > 5, 
r=1 r=1 r=1 

exist if, and only if, 0 < ab—h? + 4°(80-+-7) and of course a > 0. 


Take a = 5, so that 


(5. 
(5 


] 
] 
= 


Wr 


a | =~ | is aul 
a, = +2, a,= +1, a,=—a,— 0, 
and other solutions derived by permuting these values. The signs 
are of no importance as —b, can be written for b,. Hence 
_— | . 2.1 A2_1A2_1 A? 
h = 2b,+bs, b = b3+b3+b2+03. 
Take h = 0, then clearly 
— Bh? h2-Lh?2 
b = 563+-b5-+-0; 
has integer solutions, if and only if 5b 4 4?(80+7), ie. b is not of 
the form 4?(80+-3), and this is (5.8). 
So = (2b, —h)?+b3+-b3+03, 


> 5 S.... {5 ¢ 2 Rh2 1 Ap2 
ol 5b—h? = (5b, — 2h)? + 5b3 + 5b}, 
has integer solutions, if and only if 55—h? + 4°(80+-7). Write 

9 ~ 9 ~_9 
m = x*-+-5y?+ 52°. 
On taking residues (mod 5), this certainly has no integer solutions 
unless m= 0, +1(mod5). We can then write m= 5b—h? where 
h=0, +1, +2, and clearly we can take the sign of x so that 
x =-—2h(mod5). Hence integer solutions for x, y, z exist, if and 
only if m=0, +1 (mod 5) and m ~ 4°(80+7). This is (5.14). This 
result, which includes (5.8), can also be easily proved on noting 
m = x?+-(2y+-2)?+ (y—2z)?, 

and using (5). 

Take next a = 6, so that a, = 2, a, = 1 =a, a,= 0, 

h = 2b,+6,+5s, b = 6} +63+-02+-0?, 
whence b = 6} +b3+ (h—2b,—b,)?+0?, 
from which 
6b —h? = 2(3b, —h)?+ 3(2b, + 2b,—h)*+ 6b}. 
This is possible if and only if 0 < 6b—h? + 4°(80-+-7). Now if 
m = 2x?+ 3y?+ 62? 


* Dickson, ‘Integers represented by positive ternary quadratic forms’: 
Bull. American Math. Soc. 33 (1927), 63-70. (5.9) is proved in his paper 
‘Quaternary quadratic forms representing all integers’: American J. of Math. 
69 (1927), 39-56, pp. 42-4. 











A NEW WARING’S PROBLEM 285 


is solvable, on taking residues (mod 6), we find m = 0, 2, 3, 5 (mod 6), 
i.e. m= 6b—h? where 0<h<6. Then the sign of x can be taken 
so that x = —h (mod 3), i.e. x = 3b,—h, and as y = —h (mod 2), we 
can put y= 2b,+2b,—h. Hence the equation is possible, if and 
only if m= 0, 2, 3, 5 (mod 6) and m + 4°(80-+-7), and this is (5.15). 
In the particular case when m = 0 (mod 6), the result simplifies on 
writing 6m for m, 3x for x, 2y for y. We see then that 


m = 3x?+ 2y?+-2? 
is solvable, if and only if 6m + 4°(80+-7), i.e. m 4 2*P+1(80-+-5), and 
this is (5.4). 

The case a = 2 leads to (5.1), while a = 3 leads to a result included 
ina = 6. The results (5.2), (5.6), (5.13) cannot apparently be deduced 
in this way. 

An interesting result follows from a= 4. Then 

a,=l=a,=a,=a@ or a,=—2,a,=—a,=—a,=0. 
Suppose / is an odd number, then a, cannot be 2. Hence the equa- 


tions, (6 > 0, h odd) 4 4 
h == > b,, b = > b? 
r=1 r=1 


are possible in integers b,, if and only if 0 < 4b—h? + 4°(80+7). 
This means that b must be odd, an obviously necessary condition 
since 6,2—b, is even. Hence, if 6 > 0 and b, h are odd numbers with 
4b <h?, the equations always admit of integer solutions in b,. They 
are well known for their connexion with Cauchy’s proof of Fermat’s 
theorem on polygonal numbers, where, however, the additional im- 
portant restriction is required that b, > 0. 

Similarly, on taking h = 0, a= 7, 11, 12, 14, 15, we have the first 
five formulae in the scheme, in which the right-hand side denotes 
those numbers, and only those, which cannot be represented by the 
corresponding left-hand form: 


2x?+ 2y?+ 322+ 2ry—2yz, 4°(80-+- 1) 
2?+ 2x?— 2ry+ by, 4°(80-+-5) 
22-2224 2ry+ 2y?, 4°(80-++-5) 
22+ 3x24 2ry+ 5y*, 27° +1(80-+- 1) 
327+ 222+ 2ary+ 3y?, 4°(80-+1) 
327+ 2a?+- 2ay+ 8y?, 37°(30+-1). 
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The third is given by Dickson and the sixth is proved in §6. We 
have for the proof of the third, a = 12, 
a,=2=—4a,=43,q=0; or a,—3,a,—=1=—a;= ay. 
The first case leads to 
b = b2+-62-+.03-+4 (b,+b,)?, 


or b = 62+ 267+ 2b,b,+ 263. 


The second gives 91291591 sen 2h 1k‘ 
& b= b?-+63-+62+ (3b, +6,+6,)? 
or writing b, = 2, b, = y—2, b, = z—2, 
9 . oe .\2 | - 2 
b = «7+ (y—2x)?+ (2—2x)?+ (e+-y+z2)? 
= 4x74 2y?+ 2yz+-22?. 
Both results are included in the result that, if m > 0, 
m = x? 2y?+ 2yz+ 22? 
is impossible, if and only if 12m = 4?(80-+-7), i.e. if m = 4°(80-+5). 
6. The results of §3 are particular cases of a group of theorems. 
The essential points of the proof are 
4 
(A) the form > x? is capable of composition with itself; 
1 
(B) it represents all integers m > 0; 
(C) the number ¢(m) of representations* of m is a function of m 
which satisfies 4(m)¢(n) =A¢d(mn) where A is a constant pro- 
vided m and n are subject to certain restrictions, e.g. [m,n] = 1 
or 2andA=1. 
4 
(D) simple conditions exist for the representation of m by > 2?. 
r=2 


There are other quaternary forms having these properties. Thus 


if p, g are any constants, 
fr=ajt+pyitqitpaty (6) 
has the property of (A). For, as is well known, f = f,f, where 
= XyXQ—PY1 Yo— F122 — PN la, 
Y = XyYot Yet Gyle— he, 
% = Uy q— PY la+%%_ + Phy Yo, 
t= UyloTYy22—% Yo be. 


* Sometimes, as in my Messenger paper, ¢(m) denotes the number of repre- 
sentations of 4m; also restrictions may have to be imposed on 2,, e.g. 2, is odd. 
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It has the property (B) by a theorem of Liouville* in the seven 

cases p= 1, q=1, 2, 3 and p= 2, g= 2, 3, 4, 5. The properties (C), 
(D) are also satisfied in a number of instances. 

One illustration will suffice, say p= 1, g= 3. Then the number of 

representations} of an odd number m by the form (6) is four times 

the sum of those divisors of m prime to 3. Hence, if m, n are both 


odd and [m, n]= 1, fi)f(n) = 4f(nn). 
On considering the linear substitution 
( a,+vV(—p)yy,  Vgz+ — 
—Nqz,+V(—pg)t,, a,—V(—p)y)’ 
it is easy to prove that, in (6.1), the same solution 2, y, z, t occurs 
four times since there are only four integer solutions of 
L= 84 P4824 37, 
namely, 
€=+1,9=C€=7r=0 and »=+1,§=f=7=0. 
Since }[4f(m) ][4f(n)] = 4f(mn), 
this means that all the representations of mn by (6) are given by (6.1). 
Next, from (5.13) A= 27+ 3y?+ 32? 
is possible if, and only if, 0< A 4 3%(30+2). Hence, if a and b are 
odd and positive, 
a = a}-+-a}-+ 303+ 3a2, 
h = a,b, +a gb.+ 3ayb,+ 3a,b,, 
b = b[+-63+-3b3+ 3b] 
is possible in integers a,, b,, (r = 1,...,4), i.e. 
aX?+ 2hX Y+bY?= (a,X+b,Y)?+(a,.X+6,Y)?+ 
+3(a,X +6, Y)?+3(a,X +6, Y)?, 
if and only if 0 <A +3*(30-+2). 
The conditions a and b odd, [a,b] = 1 are not necessary. For if 
d= [a, h, b], d can be represented by (6) and the theorem is true if 


it holds for 


avo oh b 2 
gute x¥ +o Y?. 


Hence we may take [a, h, b]=1 or 2. Then the condition that 
a, b are odd can be replaced by the one that [a, h, b] A 2. 
* See Dickson, History of the Theory of Numbers, 3 (1923), 225-33, for 


various results on quaternary forms. 
3695 U 
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Hence the restriction on a, b in (6.2) is that for 8 > 1, a, h, b should 
not have 28, 28-1, 28 respectively for the highest powers of 2 dividing 
them, i.e. the form aX?+2hX Y+-bY*? must not be a multiple of an 
improperly primitive form. 
For a few numerical illustrations take h = 0, a= 5. Then 
a, = a,=—a,= 1, a,= 0, 
0 = b,+-6,+3bs, 
b = b3+-b3+ 3b2 + 303, 
whence b = (b,+-3b,)?+-63 +- 363+ 302, 
and so b = 222+ 2xy+ 8y?+- 32? 
is solvable if, and only if, 5b 4 37°(30+-2), ie. b A 3°°(30+-1). 
So h=0, a=6, a, =a,=0, a, =a,=1, 6,+6,=0 gives (5.9). 
Also h = 0, a= 2, a, =a,=1, ag=a,= 0, 6,+6,=0 gives (5.16). 
Finally a few words as to what corresponds to § 4. We may take 
some multiple of a fifth square, e.g. once, twice, etc., and investigate 
the problem as in §4. But I content myself with proposing the 
following problem. For what constant values of p,, r= 1,...,5 can 
every positive definite quadratic form (1) be expressed in the form 


aX*+2hx Y+bY? = S p,(a,X +6, Y)?, 


r=1 
where a,, b, are integers ? 











ON THE RELATION BETWEEN CORRESPONDING 
PROBLEMS IN PLANE STRESS AND IN 
GENERALIZED PLANE STRESS 


By L. N. G. FILON 
[Received 25 August 1930] 


1. In Love’s Mathematical Theory of Elasticity, § 146, is given the 
theory of two-dimensional generalized plane stress for an isotropic 
homogeneous perfectly elastic plane plate in equilibrium. Although 
the treatment there given omits consideration of body-force, it is 
easily generalized to include this case, as follows. 

The faces of the plate are taken to be z= +c, and are assumed 
to be free from traction. The actual displacements and stresses will 
be denoted throughout by w, v, w, x, yy, 22, yz, 22, Zy, and the mean 
displacements and stresses in the plane of the plate (means being 
taken with regard to z) by U, V, x, yy, Zy. As in Love’s Elasticity 
the relations between the mean displacements and stresses are found, 
on the usual assumption of generalized plane stress, that the actual 
stress 22 is zero throughout the material, to be given by: 


x = '(0U /ex+V /dy)+ 2peU /éx, (1.1) 
yy =X’ (aU /ex+eV /éy)+2ueV /éy, (1.2) 
ty = w(0U /ey+eV /ex), (1.3) 
where XN’ = 2u/(A+-2y), (2) 


\ and p» being the usual Lamé constants. 
The first two body-stress equations then give, when we take 
account of a body-force derived from a potential Q, 


04x Ox + Oty Oy +042 /02 = peQ/éx, 
ecky [0x + oxy [ey + Oy ez = peQ/oy, 


and, averaging these equations over the thickness of the plate and 
using the fact that <z = 0, yz = 0 when z = +c, we find 


Oax Ox + ety ey = peQ/ex, (3.1) 


oaty |x + dyy/dy = peQ/ey, (3.2) 
u2 
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so that the mean stresses can be derived from a stress-function y of 
x, y, in the form 
ra = pQ-+-é 2y oy”, 
yy = pQ+éy/éx2, 
ry = —O?yx/dxdy. 
If we solve back (1.1), (1.2) for 0U /ex, 6V /ey, we find that 
Eu /éx = xx—nyy, 
EeV /dy = yy —7xx, 
where E is Young’s modulus and 7 Poisson’s ratio. 
Eliminating U, V between (5.1), (5.2), and (1.3), we obtain readily 
the consistency relation 
(1+ )(0?x%x/0x®+ Oyy ey? + 20%zy /exdy) = Vi(ex+yy), (6.1) 
where V} stands for the operator é7/éx?-+-6?/dy?. 
Substituting from equations (4) into the left-hand side of (6.1), we 
obtain the equation 
(1+-)pViQ = Vi(zx+-yy), (6.: 
and, completing the substitution on the right-hand side, 

Vit Vix +p(1—7)2,} = 0, (6. 
which is the differential equation satisfied by the stress-function in 
this case. 

Note carefully that, since the plate is assumed homogeneous, the 
density p, as well as all the elastic constants, is independent of 2, y, z. 

2. It is also known (see Love’s Theory of Elasticity, §§ 145, 302) 
that in such a plate a stress-system can exist, such that, not only 
2z, but #2 and yz vanish throughout the body of the plate, the only 
non-vanishing stresses being those in the plane of the plate. Such 
a system is spoken of as one of plane stress. If, further, the displace- 
ments wu, v are even functions of z, and w an odd function of z, the 
middle surface of the plate is distorted entirely in its own plane. 
The stresses x, yy, y are then even functions of z. This may be 
referred to briefly as symmetrical plane stress. 

It is clear that symmetrical plane stress satisfies the conditions 
required for generalized plane stress. In fact, it more than satisfies 


them, since it implies the vanishing of the shears #2, yz, which is not 


necessarily required for generalized plane stress. 
Thus to every symmetrical plane stress there corresponds a 
generalized plane stress. The object of the present note is to show 
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that, subject to certain conditions to be satisfied by the potential 
of the body-forces, the converse holds good, namely, that to every 
generalized plane stress there corresponds a uniquely determined 
symmetrical plane stress, and to obtain, in a simple form, the dis- 
placements and stresses in this symmetrical plane stress directly from 
the mean displacements and mean stresses in the corresponding 
generalized plane stress. 

In §§ 15-17 some examples are given of the application of this 
procedure to certain important cases of flexure of vertically flat 
girders and cantilevers. 

3. It is to be noted, in the first place, that, if such a correspondence 
is-to exist, the potential Q of the body-forces must be independent 
of z, and not merely symmetrical in z, for no plane-stress solution 
(for equilibrium) can exist unless 0Q/éz = 0. 

This follows at once from the third body-stress equation 

0x2 [Cx + Oy2/Cy + 022/0z = peQ/éz, 
since for plane stress #2 = yz = 22 = 0. 

We may now omit the bar over the Q in § 1 and treat Q as a func- 
tion of a, y.only. 

4, From the usual stress-strain relations 

4x = d(Ou/ex+Ov/ey+ ew/éz)+ 2péu/ex, (7.1) 
Yy = A(eu/ea+ ev/dy + éw/éz) + 2pev/ey, (7.2) 
22 = N(Cu/ea+-Ov/éy + ew/ez) + 2pew/éz = 0, (7.3) 

we find easily, as in generalized plane stress, 
£x = X' (Cu/ea+-év/dy)+ 2peu/ea, (8.1) 
Yy =X‘ (Cu/éax+ ev/éy)+ 2pev/ey. (8.2) 

To these last equations we may add 

xy = p(eu/ey+év/éx). (8.3) 


We have also Nou /éa = £x—nyy, (9.1) 


Eév/y = yy—ntx, (9.2) 

Edw/éz = —y(#x+-yy). (9.3) 

All the above are identical with the corresponding equations of 
generalized plane stress. 

5. Further, the body-stress equations give (since #2 = ¥ 

0xx /Ox + Oxy /Cy = peQ/éx, (10.1) 

xy [Ox +-byy ey = peQ/ey, (10.2) 


~ 


= 22 = 0) 
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whence the stresses may be expressed in terms of a stress-function 


x, as follows: = pQ.+ 82x, /dy?, (11.1) 


Yy = pQ+-é?x, /dx?, (11.2) 
Yy = —O*y,/Oxdy. (11.3) 
These have the same form as in generalized plane stress, but now 
x, is a function, not of x, y only, but of 2, y, and z 
The equations of § 4 clearly lead to consistency relations of the 
same form as (6.1), (6.2), and (6.3), namely, 
(1+ 9) (@*tx/0x*+- 6 yy /Oy*+ 20%ky /Oxdy) = Vi («a 
(1-+-7)pVjQ = Vi(xx-+-yy), 
Vit Vixi + pQ(1—y)} = 
6. If we now otiiaae corresponding i in the two systems, 
we obtain 
tx —aXx = N'{0(u—U)/éx+-0(v—V)/6y}+-2n6(u—U)/éa, 
Gy —Iy =N{e(u—U) /ea+-2(v—V) /2y} + 2n0(v—V ) ey, 


14.1) 
14.2) 
14.3) 
14. 


~ 


( 
( 
Eé(u—U)/ex = tx—Xx—(yy—yy), ( 
Bo(v—V) joy = : YY —Yy —n(&#x— 2x), ( 
ty —xy = p{o(u—U)/ey+e(v—V)/éz}, (14 
( 
( 
( 
( 


14. P 
14.7) 
14.8) 
14.9) 


7. So far, we have not taken account of the conditions satisfied 


xx— rx— € "(xi —x)/ey?, 
yy —Yy = (x, — x) /ex?, 
ty—xy = —0*(x,—x)/exdy, 
where Vi(xi—x) = 0. 


for plane stress a=0, =0 


which lead to 
Gu/dz = —Ow/dx, dv/dz = —Ow/oy, 
and these may be written, since @U/éz = 0, eV /éz = 0, 
é(u—U)/az = —éw/éa, 
o(v—V)/éz = —Ow/dy. 
Differentiating with regard to z, we have 
e(u—U)/éz2? = 
6%(v—V)/é22 
whence, using (9.3), 
dx)(ax+Yy), 
6 /y)(xa-+ yy). 
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Differentiating (17.1) with regard to x, (17.2) with regard to y, and 
adding 

(6/6z*){e(u—U) /ex+-e(v—V)/ey} = (n/E)Vi(er+yy), (18) 
whence, comparing (18) and (12.2), we find 
(é2/é22){0(u—U) /ex-+0(v—V)/ey} = {9(1-+-9)p/ E}V3O 
= (np/2p)VjQ. (19) 
Integrate (19) twice with regard to z, and note, (i) that the term in 
z must be absent since u, v are to be even in 2, (ii) that the mean 
value of 0(u—U)/éx+-é(v—V)/éy is zero; we have 
e(u—U) |x +0(v—V) ey = (np/4u)(22@—Je*)V3Q, (20) 
or, using (14.3) and (14.4), 
(1-9) (tx + yy —ex—yy) = 4np(1+)(2?—fe*)VjQ. (21) 
Hence, from (14.6) and (14.7), 
Vi(xi—x) = {en(1+-9)/2(1— 9) }(2*— fe?) VjQ. (22) 
8. Comparison of equations (22) and (14.9) shows at once that, for 
a plane-stress solution to exist, we must have 
Vin = 0. (23) 
This, however, leads to no‘real limitation if the field of body-force 
is purely gravitational, for, since Q is strictly two- a ae 
V3?Q_ = —4yzp by Poisson’s equation, and p being constant, Vip = 
and therefore Vj{Q = 0. 
If, however, the field of force were due to a distribution of electric 
charges, (23) might lead to an effective limitation. 
The complete limitation on Q turns out, however, to be narrower 
than is implied by equation (23). 
Differentiating (14.3) twice with regard to z and using (17.1), we 
obtain 
(2? /xx®) (ee +-yy) = (0?/02%) (ex + yy — x — yy) — (1+ n)(6*/62*)(Yy—Yy), 
and, using (21) and (14.7), 
(2/0x®){n(ea+-yy)+ (1+-9)(6*/22")(x1—x)} 
= {np(1+-y)/(1—n)}ViQ. (24.1) 
Similarly, from (14.4), 
(62/Oy*){m(@x+-yy) + (1+ 9)(@/é2*)(x1—x)} 
= {np(1+m)/(1—n)}VjQ. (24.2) 
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Differentiating likewise (14.5) twice with regard to z and using 
(17.1), (17.2), and (14.8), we find, on rearranging, 
(8? /@ady){n (ex +-yy) + (1+ 9)(@/€z*)(x1—x)} = 9. (24.3) 
To discuss these equations, write for shortness 
n(%x+yy)+ (1+ )(6*/éz*)(x3—x) = R. (25) 
Then equations (24) become 
O27 R/ox* o* Rey” = {np(1+ n)/(1—n)}ViQ ) (26) 
C2 R/exoy = 0. } 
From the last of equations (26), treating for the moment 2, y as the 
only variables, we get R= F(x)-+Qy). 
The first and second of (26) then give 
F" (x) = G"(y) = {yp(1+m)/(L—n)}VjQ, 
whence clearly F’(x) and @”(y) must each be independent of x, y; 
hence V?Q is independent of x, y, and, since it is otherwise inde- 
pendent of z, V2Q. = const. (27) 
As explained previously, this condition will always be satisfied if 
Q is a gravitational potential, in virtue of Poisson’s equation. 
9. If we now substitute into (17.1) and (17.2) from (21) for x+y 
and remember (27), we have 
2.(6?/é2z?)(u—U) = {y/(1+- n)}(E éx)(ax-+-yy), 
20(0%/é22)(v—V) = {m/(1+-9)}(2/ey) (ee +99). 
Integrate these twice with regard to z, omitting odd powers of z and 
adjusting the arbitrary functions of zx, y so as to make the mean 
values of u—U, v—V zero: 
2u(u—U) = {n/2(1-+-)}(z?— 4c) (@/ex)(ax+-yy), (28.1) 
2y(v—V) = {n/2(1-+m)}(e2@—Je2\(B/ey(te+-gy). (28.2) 
These give the required corrections to the mean displacements. 
10. To find the stresses, substitute into equations (14). The pro- 
cedure is simplified if'we write (14.3), (14.4) in the form 
rx —Xx = 2ne(u—U)/ea+{y/(1+)(ex+yy—Xxe—yy), 


yy —Yy = 2pe(v—V) /ey+{n/(1+y)}(ex+yy—®e—yy), 
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and now substitute on the right-hand sides from (28.1), (28.2), and 
(21). The result gives 
tx — HK = {n®p/2(1—n)}(z2?— 4c?) VIQO+ 

+{n/2(1 + n)}(2?— fe*)(0?/6x*)(wx+-yy), (29.1) 
yy —Yy = {n°p/2(1—)}(22@—4e?)VjQ+ 
+-{n]2(1-+m)\(22—4e2)(e2/ey2(e+-yy). (29.2) 
For the shear <y substitute from (28.1), (28.2) into (14.5) direct 
and we have 
ty —xy = {n/2(1+-m)}(2*—}e*)(0?/Oxey)(ex+yy). (29.3) 
These give the required corrections to the mean stresses. 
11. To find the transverse displacement w, we substitute from 
(28.1) and (28.2) into (15.1) and (15.2). 
— 2pdw/ex = {nz/(1+-)}(e/ex)(ex+-yy), 
— 2pow/ey = {nz/(1+-y)}(6/ey)(Fx+-yy), 
— —2uw = {yz/(1+-y)}(ex+-yy)+H(2). 
To determine H(z) we use equation (9.3). Thus 
n(tx+-yy) = n(Xe+yy)+(1+-y)H'(2), 
which, from (21), leads to 
{y?p(1 + 9)/2(1—m) }(22— $e?) VQ = (1+) H(z). 
Integrating this, and remembering that w, and therefore H(z), must 
be odd in z, we have 
H(z) = {y2p/6(1— n)}(z8—c?z) ViQ, 
so that 

2uw = —{n*p/6(1—n)}2(2*@—c?) ViQ—{nz/(1+-9)}(ex+yy), (30) 

which gives the third displacement and completes the solution. 
12. It is desirable to obtain explicitly the function x,. Referring 
back to equations (26), and having regard to (27), they clearly lead to 
R= }(a®+y?){np(1+ y)/(1—n)}VjQ+-arbitrary function of z. (31) 
(25) then gives us 
(1-+-)(@/@2?)(x1—x) = R— ae Seige. 
= a*+-y*){ne(1 + 9)/(1— 9) VjQ— (ex + yy) + 
+-arbitrary function of z, 


the middle term on the right-hand side having been absorbed into 
the arbitrary function, since (21) and (27) together show that the 
term in question is a function of z only. 
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Integrating now the equation last written, and using the condition 
that y,—x is to be an even function of z whose mean value is zero, 


we obtain 
X1—X = H(2e*+-y*){ne/(1— 9) }(P— 3?) VjQ— 

—{n/2(1+ n)}(P—fe*)\(ex+yy), (32) 
the arbitrary function of z being finally omitted, since it clearly 
would contribute nothing to the stresses and may therefore be 
dropped without affecting the solution in any way. 

It is easily verified that, if we substitute from (32) into (14.6), 
(14.7), and (14.8), bearing in mind (6.2), these equations are satisfied. 
The solution is therefore consistent, since all the equations (14) have 
been shown to be satisfied, and it is clearly uniquely determinate. 

13. We have therefore shown that, provided the force-potential 
Q is such that VjQ = const., there exists always a symmetrical plane 
stress corresponding to a given generalized plane stress, and that the 
actual displacements and stresses in this symmetrical plane stress 
are derived directly from the mean displacements and mean stresses 
by adding corrective terms, as given by equations (28) and (29); and 
these corrective terms are known when the single quantity zx+-yy 
has been determined in’ the generalized plane stress, and do not 
depend directly on the boundary conditions. 

Moreover, if we refer to (30), we see that 

24(W4p—We) = —{2ne¢/(1+-m)\ te +9y), 
so that «x+yy = (E/2nc) x (lateral contraction of the plate at the 
given point). 

In photo-elastic investigations of stresses in plates considerable use 
is made (e.g. by Professor Coker) of such measurements of contrac- 
tion or expansion by means of a lateral extensometer. The data given 
by such measurements, provided a formula can be fitted to them 
reasonably well, so as to enable the differentiations in equations (29) 
to be carried out, will enable the experimenter to estimate the de- 
parture of the actual stress in the plate from the average, a'point of 
some importance. 

14. The most frequent and useful case occurs when Q is the 
gravitational potential of masses at a distance, usually the earth. 
In this case VjQ = 0, and, if we write 


{n/2(1+-y)}(P—fe*)(ex+-yy) = $, (33) 
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the solution assumes the extremely simple form 
x1=x—-¢, (34.1) 
2u(u—U, v—V, —w) = (0/@a, 6/ay, @/éz)d, (34.2) 
(tx — Fx, Ly —Fy, yy —Jy) = (02 /ex®, 2/exdy, 2/2y2)d. (34.3) 
15. We will now take a few examples of the applications of these 
results. 
Consider the generalized plane stress given by 
x = My /6I, Q=0, 
I being the moment of inertia of the rectangle y= +b, z= +e 
about the axis of z. 
This corresponds to 
tx=My/I, yy=xy=0, 
and therefore to a flat girder of thickness 2c and height 26, bent in 
its own plane under a uniform bending-moment M. 
It is then found easily by the usual methods that 
U = Mexy/EI, V = —(M/2E1)(x?+ ny’). 
We have now 
$ = {n/2(1+-9)}(?—}e*)(My/Z), 
whence 
ex = My/I, 
0, 


~ 


yy = yy = 9, 
u = U+(1/2u)0¢/ex = U = May/ET, 
v= V+(1/2u)ob/ey = —(M/2B1)fa2+ nly*—2?+ Je}, 
w = —(1/2p)ép/ez = —(M/E1)yz. 
These are the well-known displacements of Euler-Bernoulli flexure, 
and it is at once seen that, in this case, they apply to a beam of any 
cross-section, since the stress across any element parallel to the axis 
of x vanishes. 


16. Take next the case where 


x = (W/21)(4ay?—b22y), Q=0, 
the symbols having the same meaning as in the last section. 
This leads to 
te=Way/l, yy=0, xsy=(W/21)(b?—y’), 
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corresponding to a flat cantilever, whose plane is vertical, and which 


carries a total transverse load W at «= 0. 

Here $ = {n]2(1+-9)\(e2—3e2)Way/L, 

and Le = £2, yy = YY; 

iy = Fy +{n/2(1-+ n)}@2—Je2)W/I - (35) 

= (W/21)[b2—y?+-{n/(1+-n)}(22—4e2)]. |) 
In this case the top and bottom faces are under a residual traction 
ty = {Wy/21(1+y)}(2?—4c*), 

which is, however, self-equilibrating at every cross-section, so that, 
applying de Saint-Venant’s principle, the stresses inside the body of 
the plate should be given to a sufficient approximation by (35). If, 
in fact, the well-known solution for flexure with shear of such a beam 
of rectangular cross-section be examined, and terms which become 
exponentially small at a distance from the top and bottom be 
neglected, the approximate stresses deduced are found to be identical 
with those given by (35). 

The displacements in this case also are readily calculated as before, 
if required. The only laborious part is the preliminary calculation 
of U and V, but the algebra is quite straightforward. 

17. Finally, to illustrate the case where there is a body-force, we 
will consider a problem for which 

x = Aa®y+ Ba®ye+ Cy>+ Dy’, Q = gy. 
Since Viy = 0 (for V?Q = 0), we must have C = — B/5, so that 
x = Axty+ Bx*y®— hy’) + Dy’. 
Choose now A and B so that the mean stresses #y, yy both vanish 
over y= +b. This requires 
0*y/0x? = —pgy and 0*y/dxéy = 0 
when y = +b, leading to: 
A+ Bb? = —h0q9, A+3Bb? = 0, 
or A = —$pg, B= }(pg/b?), so that 
x = 4F(pg/b?)(a?y— 3x2yb?— Ly5) + Dy’. 
This gives the generalized plane stress for a heavy uniform girder of 
density p, height 2b, and thickness 2c, under its own weight. 
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The mean stresses are given by 

ir = (pgy|2b*)(2b2-+ 3x2 2y2) + 6Dy, 

Jy = (pgy/2b*)(y*—b°), 

&Yy = (3pgx/2b*)(b>—y?). 
To find the total transverse shear S per unit thickness and the 

bending-moment M per unit thickness, we have 
b 
S= [ ty dy = 2pgbar, 


—b 


b 
M = | rx y dy => — pgb(x*+- +b?) —4Db?. 
—b 


If the girder is simply supported at its ends, which are distant 
a from the plane x = 0 of symmetry, then 
4Db° = —pgb(a*+- 4b"), 
and M = pgb(a?—zx?). 
The stress-function then becomes 
x= (pg/40*)(x%y®—a*y?— 32tyh?— hy spy) 
and x = (pgy/2b?){3(x2—a*) — 2y?+- 9b}. 
Hence rx +-yy = (pgy/2b?){3(x? —a?) —y?+ 1b}, 
and = {npgy/46%(1-+-9)\(2*—4e*{3(a*—a*) —y2 + $09, 
so that, in the corresponding plane stress, 
ta = xx +{3npgy/2b>(1+ n)}(2?—}e*) 
= (pgy/2b*)[3(a°—a*) —2y?+ $b? + {3n/(1+-9)}(2?—$e*)], 
Jy = Yy—{B3npgy/2b*(1+ )}(2?—fe*) 
= (pgy/2b*)[y°—b*—{3n/(1+- )}(2*—$e*)], 
ty = fy +{3npgx/2h*(1+ y)}(z?— he?) 
= (3pgx/2b*)[b?>—y?+-{n/(1+-9)}(°—}e*)]. 
These should give close approximations to the actual stress at any 
point in such a girder at a distance from the top and bottom faces. 
In this case, however, no exact solution for a heavy rectangular beam 
is available for comparison. 
The displacements, in this case also, are not difficult to obtain in 
finite terms, but they lead to long algebraic expressions and it appears 
hardly worth while to write them down. 








/' THE EXPONENTIAL FUNCTION IN 
LINEAR ALGEBRAS 
By P. DIENES 

[Received 16 July 1930; in revised form 28 August 1930] 
1. THE introduction of complex numbers was chiefly suggested by 
the problem of determining the zeros of polynomials. Some integral 
functions, such as e*, have no zero in the field of complex numbers. 
This fact suggests the following question. Can we generalize the idea 
of number to such an extent that the exponential function may have 
a zero in the extended field ? 

We shall prove in this Note that the exponential function has no 
zero in the linear associative algebras to a finite base, and that it has no 
self-associative zero in finite non-associative linear algebras. This 
result extends to a large class of algebras to an infinite base. In 
particular, the exponential function has no zero in the tensor algebra 
of relativity theory and it misses only singular tensor values which 
do not divide some tensor. Moreover, in Hilbert’s* algebra of infinite 
bounded matrices, so important in atom mechanics, the exponential 
function has no absolutely bounded matrix zero. 

Since we have to define the symbol 

1+ 

7 

for numbers x of a general linear algebra, it is necessary to suppose 

that 1 is a number of the algebra, i.e. that the algebra has a modulus 

é, such that if x is any number of the algebra, xe) = egx = x. Ob- 
viously the algebra must also contain a zero element. 


na 
> x"/n! 
1 


v 


2. Construction of an algebra to an infinite base 
The numbers of an algebra to the base é, = 1, @,,..., €,,... are, by 
ie) 
definition, the forms x = > 2,e;, where the components 2; are real or 
i=0 
complex numbers. Addition is defined as addition of the correspond- 
ing components. Multiplication is supposed to be distributive with 
respect to addition and is defined (a) by the multiplication table 


io 2) 


a 
e;e; ~ -) Mijas (1) 


* D. Hilbert, Grundziige einer allgemeinen Theorie der linearen Integral- 
gleichungen, Leipzig (1912), pp. 128-9. 
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where m,;,, are arbitrarily given real or complex numbers, i.e. by 
i a) ao 
ry= > (3 m1] Cs (2) 
k=0\i,j=0 
and (6) by the convention that if c is any real or complex number, 
Ce == B6= 5 Cxg,. (3) 
i=0 
Zero is the number of the algebra whose components are all zero. 
We see that 02 = 20 = 0. 
We also notice that multiplication does not lead to a number, 
unless the infinity of double series 


; 2, a; YM 45, (k = 0, 1,...) 
i,jJ= 


are all convergent, so that, in every algebra to an infinite base, there 
are numbers whose product is meaningless, unless, for every k and 
for i>N, 7 >N, we have m,;,—0. In fact, we can choose 2;, y; 
such that x; y,;m,;, = 1 for an infinity of suffixes i and j. 

The necessary and sufficient condition that the multiplication of 
units should be associative is that (e,e;)e,, = e;(¢,e;.), i.e. 


i 3) wo 
> MijNoaK = > Mod jhe (4) 
a=0 a=0 


for every i, j, k,l. In an algebra to an infinite base these conditions 
are not sufficient to make multiplication associative in the whole 
algebra. In fact, the detailed form of the equation (xy)z = 2(yz) is 


oe U5 YjeKRM ij MN ak = >> 4% Yi=RM jo Niod> (5) 
ak ij ai jk 


and unless the order of summation can be inverted, (5) is not neces- 
sarily a consequence of (4). 

For example, in the infinite matrix algebra to a double array of 
units e;; with the multiplication table 

©57ea = * = ¥ ; (6) 
multiplication of units is associative, but for suitably chosen coeffi- 
cients, the equation 
(> Gijig X Opes) X Cmn&mn = DX Vijeiz (X Oxelnr X Cmn&mn) 
i.e. z (> Viabap)“Bn = > Gig > bapepn 
| 0% a B 

is not satisfied. This remark seems to be necessary, for in papers on 
quantum theory the associative property is taken for granted also 
for non-bounded infinite matrices. 
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3. Power series 


2 
We say that a power series > a,x”, where the coefficients a, and 

n- 0 
the variable x belong to a given algebra, is convergent if all the 
products a,x" => (a,x”"),e; have a sense and if all the ordinary 


t 
infinite series > (a,2”);, (¢=0,1,...00), are convergent. We readily 
a 


see that, if c is any real or complex number, 
[a,,(cx)"}; - c"(a,x"); 
as soon as (a,,2”); exists. Hence: 


THEOREM I. Jf > a," its convergent, > a,(cx)" is absolutely con- 
vergent provided the absolute value of the complex number c is less than 1. 
We shall now examine the Cauchy product of two infinite series. 
If the product ~ : . ; : 
¥ ~* S 
7a ay a bn) moe MN ik a (@,,)i > (b,); 
n n yk ij n n 
exists, we have 
>a.> +. } l 
(> a, >. bn) = > Mix > [(a9);(6,,);4 (41) (O,-1)j; +--+ (Gp)e(o),], 
n n Tk ij n 
provided the rearrangement is legitimate, e.g. when the two series 
> (@,,); and > (6,,); are absolutely convergent (or at least one of them). 


n 7” 


If also the order of the two successive summations, viz. with respect 
to n, and with respect to 7, 7, can be inverted, we obtain 


= y Stee a8 ih. | 7 
(2 Ay a Pn) meer [ab, I a,b, 1 | +e FOndoly 
nt q n 


vt 


5 oh. oe sy a 
> 4, > 6, => [aob,, +4,5,-1 +... +459] 
nt vt n 


which corresponds to a well-known theorem on infinite series with 
real or complex terms. 

The last inversion is obviously allowed if i and j assume only a 
finite number m of integral values, i.e. the base of the algebra is finite. 
Hence: 

THeoreM II. If, in an algebra to a finite base, both > a, and > b,, 
are absolutely convergent, the Cauchy product of the two series is abso- 
lutely convergent and we have 


P a > b, — > (ab +4ayb,,. 1tb++s + Ando). (5) 


naw n | 
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4. The exponential function in finite linear algebras 
We say that 2 is self-associative if 
a = zo? (1) 
whenever n’-+-p’ = n+p. Take a self-associative number z of a finite 
algebra and consider the series 


= x" mn) 9 
> a -> 1D (x es (2) 
n=0 
Denoting by M the largest of the numbers |m,;| and by p the largest 
* Oe \e");| <plm*p My", (3) 


which establishes the absolute convergence of all the m series of the 
right-hand side of (2) and thus proves 

THeEoReEM I. In a finite linear algebra the right-hand side of 

E(x) =1+)> x"/n! (4) 

defines the exponential function for every self- Bree number x of 
the algebra. In a finite associative algebra, E(x) is defined by (4) over 
the whole algebra. 

The exponential function so defined is continuous if we define 


continuity by the condition that f(«)—/f(x) for every sequence 
x”->« fie. lim(a™); = x; for every i]. In fact, since (2) consists of 


no 


a finite number of infinite series, uniform convergence in any finite 
domain of the algebra is established by (3) and, as usual, continuity 
follows from uniform convergence. 
The direct result of applying Theorem II, § 3, to the exponential 
function is 
E(rx) E (Fx) = E(rx+F2), r,? real or complex, (5) 
provided z is self-associative. A particular case of (5) is 
E(—2x) E(x) = E(0) = 1. (6) 
But #(—2) is an ordinary (non-infinite number) of our algebra. There- 
fore, if we suppose that E(x) = 0, (6) leads to an absurdity. Finally 
« is certainly self-associative in associative algebras. Thus we have 
THEoreM II. The exponential function has no zero in finite linear 
associative algebras. The exponential function has no self-associative 


zero in finite linear algebras. 
The question arises, ‘Does the exponential function miss values 
3695 x 
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other than zero?’ Theorem II, §3, applied to > #”/n! and > #"/n!, where 

x and & are commutable, i.e. x% = x, easily leads to the equations 

E(x) E(#) = E(x+2), if 7% = xz, (7) 

E(nx) E (px) = E(rx) = E(n’x)E(p'x), ifnt+-p=n'+p’=r. (8) 

It follows that if x is self-associative, then also E(x) is self-associative. 

Also, if E(x) =a, then a” = E(pz), ie. if a is a value of E(x), then 

all the positive powers of a are assumed by E(x). Hence: 

THEOREM III. All the values of E(x) assumed for self-associative x 
are self-associative. The exponential function misses all the nil-potent 
numbers of every finite linear associative algebra. The exponential func- 
tion misses all the self-associative nil-potent numbers of every finite 
linear algebra. 

The equation 

: iz, 2 x ,(x\ |2 
E(x — ce “T E(-) 
q q q L Nee 
leads to the definition of gth root. Ifa is a value of E(x), a= E(z) 
say, we put e 1) on = 

: | a4 — K(#/q), ava — E(px/q), (9) 

which, by the continuity of E(x), leads to the definition of a real (or 


complex) power r of a, 
a’ = E(rz). (10) 


The power so defined possesses the two fundamental properties of 


indices, viz. q’a® = art, (a’)* =a’, (11) 


If b = E(g) and &, 7 are commutable, in which case a and b are them- 
selves commutable, we also have 
a’b" = (ab)’. (12) 

The extension of the results of this article to algebras to an infinite 
base is a delicate problem. We shall discuss the exponential function 
in two different types of infinite algebras used in relativity theory 
and atom mechanics respectively. 

5. The exponential function in tensor algebra 

Consider an algebra to an infinite base where the multiplication 

table is replaced by the convention that all the products 
€;e;, (€,€;)e,, = €;(€;€;,) = e€;€,, ete. 


of the given set of ‘fundamental units’ ¢,, é,..., é,, are fresh units 
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of the algebra. Thus, besides ay + S a,e;, the algebra contains num- 
i=1 


bers of the form 


m m 


= iseie ByOgyee 2 Miamtslige Piper 
< 


tke 
where e,e; and e,e; are different iol. In this section we shall make 
use of the summation convention, i.e. we shall suppress > if the 


same suffix occurs twice in a product. Thus s Misi will be denoted 
i, j= 


by a,,e,e;. The general form of a number of this algebra is an infinite 


polynomial in é,,...,e,, with real or complex coefficients: 
A = Ag+ pl; +Ojj00j;+--- Oj, 4:01, Cit 
B= bot bye, +b j6:¢;+- +05, iC: Cag te 
Sum and product are defined by putting 
A+ B= dy+bo+(a;+5,)e;+ (a,;+5;;)ee;+...+ 
+ (44,0, +5, tig te (1) 
AB = agbyt (gb; +4 jbo )e;+ (gbj; +b; +; jbo ee; +--+ 
+ (gD 5-4, +0405, tg FO, Doli Cyt (2) 

We notice* that the product always exists, that 0A = AO0=0, and 
that if AB = 0 then either A or B (or both) are 0. 

The algebra so constructed will be called the tensor algebra to the 
fundamental base ¢,,...,¢,,, for the tensor algebra of relativity theory 
is of this type, inaly restricted, however, to homogeneous forms. 

TueoreM I. If, in a linear associative algebra, xy = yx and n is a 
positive integer, then 

(wy) = 2" cta" ly +... pepentyk tty (3) 

Proof. Multiplying xy = yx by x and y on the am and left we see 
that ay = xyx = yx®, xy? = y*x. In general a”y* = ykx™. Therefore 
any proof of the binomial theorem for real or complex x and y 
extends without the slightest modification to any pair of com- 
mutable numbers x and y of a linear associative algebra. 

THEOREM IT. The exponential function 

E(x) = > x"/n! (4) 
has a definite tensor value for every tensor x. 


* See P. Dienes, ‘Sur la structure mathématique du Calcul Tensoriel’: 
Journal de Math. Ser. 9. vol. 3 (1924), pp. 79-106. 


x2 
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Proof. Putting w=—a,4+X X = x¢;+2 


A, Cj... 
we have, by (3), 


ij 


n(n— 
Nn yn | m-1Y¥_| 
z* = = Ly TNX X + = 


ee ae 
)ym-2N24 4 Xm, 


for 2) is commutable with every tensor. 
Putting xX? = Fies.-t...+X? 


- D  o 
aa. eZ Po I 


X* == Zé 


406,008 bee 


‘ | 
we see that the coefficient of e; ...e;, in #"/n! is 
nl n-k Xk 
; XO r | . cool | xO : Xi ty (6) 
We ta are en a es 
(n—1)! ‘t""* " (n—2)! 2! (n—k)! ki ’ 
i.e. it consists of k terms (when n > k), the summation with respect 
to n affecting only 2»). Therefore, when we sum with respect to n, 
the coefficient of e; ...e;, in > x"/n! will be the finite sum 
€: rk 


Loy ao igetk rLy ~~ 144m —_ § BY»)? ad 
CPL: 4 emo — 3 +... em il "= {E(z)}; in (7) 


which proves the theorem. 

In particular, the terms without ¢,,...,e,, in > #”"/n! lead to e” as 
the Oth coordinate of E(x). But (x)= 0 means that every coordi- 
nate of E(x) is zero. Since the 0th coordinate cannot be zero, we have 

Turorem III. The exponential function has no tensor zero. 

We also notice that 

E(x+y) = E(x) Ey), if ry = ye. 
Left-hand division of A by B requires the solution of the equation 
i.e. the equations ee 
oo 0 
bit thor; = a; 


ij 


. “ | a _- 
Dj ;%y+b;2;+-bor;; = a 


Ds ttot Op, tg Vig t - +b o%;.,...é; = i, iy 
Therefore, if b, 4 0, left-hand (and right-hand) division by B is pos- 
sible and the result is unique. If 6, = 0 and a, 0, the division of 
A by B is impossible. 
We see thus that if b,=0, there are tensors A such that the 
division of A by B is not possible (B is a singular tensor), and if 
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b, #0, division of A by B is possible for every A (B is a regular 
tensor). The exponential function assumes only regular tensor values. 

On the other hand, it assumes all such values. In fact equating 
(7) to any given a,,_;,, the values 2, %;,%;;,...,%;,_4, Satisfying those 
equations are successively determined by the equations them- 
selves since X7;,... X_;, contains only 2;,2%;;,..., and the 
coefficient of a; _;, is e* 40. Therefore 


iy..0p-y 


THEOREM IV. The exponential function assumes every regular tensor 
and misses all singular tensors. 


6. The exponential function in matrix algebra 


Matrix algebra is specified by the multiplication table (2.6). We 
y, with Hilbert, that the real numbers x; form a wnit vector if 


sa 
ea 
~~ 


x? = 1 and that the matrix a,,e;; is bounded if 


0 and 
| S a,;0,x; <M for every unit vector. (1) 
ij=0 
The least number // satisfying (1) is the bound of the matrix. 
Hilbert proved (a) that the product matrix a,jb;,e,, is bounded by 
the product MN of their respective bounds; (b) that multiplication 
of bounded matrices is associative. 
It follows from (1) that 
la iz\ <M, (2) 
i.e. the components of a matrix are bounded by the bound of the 
matrix (but the converse is not true). The two properties of multi- 
plication established by Hilbert show that every power x” of any 
bounded matrix x bounded by M is a perfectly determinate bounded 


matrix bounded by M”, i.e. by (2), 
I(2" ap] <M”. (3) 
We thus see that the series 1+ > a"/n!, where 1 is the unit matrix, 
converges for every bounded matrix x. Therefore 
Turorem I. T'he exponential function is defined for the whole algebra 
of bownded matrices by 
E(x) = 1+) 2"/n! (4) 
and all its values are bounded matrices. 
Now consider two series 


> a, and > b, 
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of absolutely bounded matrices, i.e. such that A;;, = |(a,),;| and 
Bi jn = \(b,)ij{ are also bounded matrices bounded by A, and B,, 
respectively. If we suppose that both series (5) are absolutely con- 
vergent, we have 


(> 4, > On) ix ee 


> (@, ij > (b,.) sx 


nm 


Ze z [(49)ij(On) jut (41) ij, -1)jn ++ + (@y)iz(Oo) 5x1- 


nm 


S 
j 
~ 


Since, by hypothesis, a series remains convergent when we replace 


a, and b, by A,, and B,, the series converges absolutely, and thus, 


U n? 


inverting the order of summation, we obtain 
> ay > b, _ > [ayb,, +a,b,, 1 Pees +4,bo]. (6) 
n 

Since bounded matrices are associative, they satisfy (4.1), and thus, 
by the repetition of the reasoning used in the case of finite algebras, 
we obtain H(—wx) E(x) = 1 leading to 

THEOREM II. The exponential function has no absolutely bounded 
matrix zero. 


7. The exponential function in bounded algebras 

The main difficulty in examining the exponential function by our 
method in infinite algebras consists in the inversion of the order of 
summation on the right-hand side of (3.2). We can justify this inver- 
sion by a convenient restriction of the multiplication parameters 
m;;, or by restricting the series >a,, > ,, i.e., in the case of power 
series, the values of x, or, again, we might make use of all these 
restrictions. 

As an example we shall consider ‘bounded’ algebras, i.e. such that 
>, |mijx| = My, <M, (1) 


a 


t,J 


and we shall restrict x and the coefficients a,, b,, to ‘bounded num- 


n? 

bers’ of the algebra in question, viz. such that 
og | | koe 

jet | P> \(4,.)i\ < A,, \(6,,);| —_ B 


vi 


It follows that \(a,x");| << A, p"M", 


— 1 
so that, if lim |VA,,| = —, 
o that, 1 ce i * R 


the series > a,” will be absolutely convergent if Mp < R. The value 
R/M is a kind of radius of convergence for >a,2". If R’ is the 


n* 
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corresponding number for > 6,2" and R is the lesser of the two, the 


series > pA, B,+A,B, 1+... +A,B,) = C(p) (5) 


converges if p< R. 
On the other hand, if x is self-associative, commutable with all the 
coefficients a, and b,,, and if |a;| <p, we have 


( z a,x" > b,x” deel 
n 


= |> Mizz > (Go) (Oy, 2"); + (GX) (b, 12"); +--+ (Gy, 2") (09) 5] 


ij n 
< > Ms j4.| > (ApoB, + A,B, 1+... +A,Bo)p” < MC( ), 
ij 


by (5) and (1). The triple series being absolutely convergent, we can 
rearrange it at will. Thus we obtain 

THEeoreEM I. If x is a self-associative number of a bounded algebra, 
commutable with ali the coefficients a, and b,, and if R and R’ are the 
radii of convergence of Sa,x" and > b,x" respectively, and R is the 
lesser of the two, the Cauchy product > (agb,,+-a,b,,_4+...+@,b9)a” exists 
for every x such that \x;| <p <R and is equal to Sa,2" >b,2". The 
product so obtained is a bounded number. 

The direct result of applying Theorem I to the exponential series, 
where the coefficients are positive and FR infinite, is 

TueoreM II. The exponential function is defined for every self- 
associative bounded number of a bounded algebra, all its values so defined 
are self-associative bounded numbers and it misses all the self-associative 
nil-potent bounded numbers of the algebra. 

The last part of the theorem follows, as usual, from 
E(rx) E(Fx) = E(rx+Fx), showing that with a every power of a is 
assumed by E(x). 

For associative bounded algebras the last result reduces to the 
simpler form: 

TueorEM III. Jn associative bounded algebras the exponential func- 
tion exists for every bounded number, and all these values of E(x) are 
bounded numbers. E(x) has no bounded zero and misses all the nil- 
potent bounded numbers of the algebra. 
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1. Iv has been proved by Landau* that, if f(z) is a function which 
is analytic throughout the interior of the circle |z| = 1 and which is 
expansible in the Taylor series 


f(z) = ag +a,2+-a,27+..., 


and if |f(z)| <1 whenever |z| <1, then 


\@9-+€,+45-+...+a,|<G,, 
where be \\* 1.3)\2 ‘1.3.5...(2n—1)\?2 
i (5 (x a (sa8 sania 2n 


Landau has also shown that, corresponding to any given value of 
n, a function can be constructed such that |ay+-a,+a,+-...+a,| is 
actually equal to G.,. 

It is evident from Cauchy’s theorem that 


: l 7 oe 1 ) 
Ay +a,+a,+...pa, = oi EGE tat3t-+= aj dz, 
* 


| 
~ ~ 


where the contour C is the circle |z| =r < 1. 


Hence eo 
; l e : p, rh 1p(n +1)i0 ] ; 
Ayta,t+a,+...+a,|=|— | fre) a dé 
. - ’ ¢ Qa * rnentO (rei? — 1) 
p l ° ; ; rr t1e(n 1)i¢_ ] 
<— | /f(re'®) a dé 
2a | rren? (re? — 1) 


] e yr 1p(n Hid__] 


ren O(ret? nicl ) 





dé. 
5) 





— 


* Archiv der Math. und Phys. (3) 21 (1913), 42-50, 250-5; Ergebnisse der 
Funktionentheorie (1929), 26-9. 
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Since this inequality holds for all positive values of r less than 1, 


we must have 
prtig(n +1)i0__ 1 


== 
(M71 a,+a,+.. -+a,|< fee a Qn = meni (rei — 1) 


eee 
Pas. i sin }(n+1)0 do 


2a sin 30 
oe f /sin(n-+1)0 de 
7 sind | 
—ho 
2 fm (n+1 6 39 — 
sin 0 oy 


where p, denotes Lebesgue’s constant* defined by the equation 


2 | 1 dé. 


sin 6 





P, = — 
7 
0 
When we give dp, 4, @, ..., @, the special values whose sum is 


numerically equal to G,,, it is evident from these inequalities that} 
G, < Pin: (1) 
By taking n = 1 in this inequality, we infer that 7 < 3-2 
Analysis of the type which has just been given shows that, if F(@) 
is any integrable function such that |F(@)| <1, then 


| Hj i 
1 e(n+1)i__] . : 
i | FO Saea1) < Pin 
7 





and it is fairly evident that, corresponding to any given value of n, 
~ the bound p,, can be attained by choosing the function F(@) in an 
appropriate manner. 

The last integral, however, has an upper bound G,, which is less 
than p,, when F(@) is the limit, as r>1, of an analytic function 
f(re®®). It is consequently of some little interest to investigate the 
relative magnitudes of G, and p,, as an estimate of the effect of 
prescribing that F(@) should be the limit of an analytic function. 

* Lecons sur les séries trigonométriques (1906), 86-7. In work on Fourier 
series it is adequate to consider p, where n is an integer; but the properties 


of py, Which will be required here are similar to those of p,. 
+ So far as I know, the inequality (1) has not previously been noticed. 
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It is evident that G = py) = 1; and, when 7 is large, 


1 4 
G,, ~—logn, Pin = — logn+c-+o(1), 
7 7” 


where c is independent of n; these formulae are due to Landau and 
Fejér* respectively. 

The results which have just been stated indicate the possibility of 
the truth of the inequality 


! 4 
Gn < Pan < -G@. (q:= 0, 1.2.5; 
7 


part of which has already been proved, and further that, when n 
assumes integral values, p,,,/G,, may be an increasing function of n, 


so that 4 
Pin — Pyin+t) — * (n = 0, 1, 2,...). (3) 


Y bins ul 
G G, +1 w 


n 


€ 


2. The conjectures (2) and (3) prove to be correct, and they will 
be established in §§ 5, 6; before investigating them, it is convenient 
to make a digression by examining the asymptotic expansions of p,,, 
and G,, which, I think, have not been published previously. 

Take Szegé’s formula} for Lebesgue’s constant 


ws 1 f8,1.,1 l 
== oe Neaft Sande a a" 4 
Pin pa mao -— a | 2m(n eye |g (4) 


if (x) denotes the logarithmic derivate of the gamma function, 


we have 


7 


oa 
; P 1 
: f ae i 1)\% 
1 _ ————_—— 5Us( 27 m--s3 5) 
Pin 2 4m? 1 {ah( T >) b(3)} 
m 


4 << < (mn+m-+4) 
—-=4(3)+ 5 > ; =. 
7 1° 4m*—1 
m=1 


Now, in the asymptotic. expansion 


=. (—y-AB / 1 1 
mn+m-+4)~ log(mn+m)- Zz (—) r(1- ae 
H( 2) Somer, a 2r 22-1] m2(n +1)?" 


the error due to stopping at any term is of the same sign as, and is 


* Journal fiir Math. 138 (1910), 22-53. 

+ Math. Zeitschrift, 9 (1921), 163-6. It is worth mentioning that this was 
obtained very ingeniously by expanding |sin d| in a Fourier cosine-series, 
deducing that a) oe 

Isin | = 8 Ss sin?md 
T m= 


writing (n+1)@ for ¢ and integrating term by term. 


1 4m? —1’ 
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numerically less than, the first term neglected; and, since (4m?—1)-* 
is positive, it is easy to see that when we substitute the asymptotic 
expansion of %(mn-+m-- }) in (4), we obtain the asymptotic expan- 
sion of p,,, with the same properties for the error. That is to say, 

4 ; logm . 
Pin ~ =| login+1)+ {2 mtn ¥4)| + 


q 


+ > Sa ze) > a rm} arp 


r=1 m 





so that 


4 y- - 
— (n+1) +g +e ped 1 =| (5) 


io 2) 
~ logm 


4m2—1 
m=1 


c= FH(1— e =) > , (r =1,2,3,...). (7) 


4m?—1)m> 


where (3), 


The constants c, (r > 1) are ees in finite terms; thus 


i) wo 


1 1 
mip = 2 amt ony —~ Gai Gm} 


“pes > el 
2 


so that : 
92r-1_.])B em - 
oul el > Ia | 2,3,..). (7a) 
m=1 \~ : 
The constant c, is expressible in a form suitable for computation 
as follows: 


logm _ > logm _ ve) 
4 <3 (sre 2m--1 


=~ (log(m+1) logm =>; 1 ” (2m+1)+1 
~2o areg > Im+1 8 Qm+1)—1 


ss ieee 
= sid erie") ea 


2) 














1 
Ww here U, 2r+2 — -> (2m+ 1)?r+2 . 


m=0 
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The values of U,,, 


U,,..—1 decreases sufficiently rapidly, as 7 increases, to make the 


io 
, = Usse—1 
computation of 2 > 2 3 i an easy matter. 
. r a 


2 have been computed by Glaisher,* and 


r=U 
We thus get 
"9 = 2-44132 38136 94835, *4 = 0-00275 82678 40582, 
c, = 0-02958 88277 62644, | c; = 0-00504 16328 46519, 
Cy = 0-:00492 82989 42233, *, = 0-01405 56865 05171, 
Cs = 0-00257 11773 16860, | c, = 0-05554 94529 62768. 
I have calculated p,,, when n = 9 by means of this asymptotic expan- 
sion and also by direct calculation from Lebesgue’s integral, using 
Andoyer’s fifteen-figure tables of the circular functions. The results 
agreed to fifteen significant figures, giving 
Pog = 1°92 35823 652 
The asymptotic et of io is a little more troublesome to 


investigate. Take Ramanujan’s formulay 


GC. - (I'(n+8))?( 1 -(5) l Fi 
"~~ ([Pin+1)} (n+1 5 (2) n+2 © \2.4) n+ 


whence it follows that 


| (4 -(5 ) - +(54) #+-| ke 


1 
: Kk2+) dk == * | K'(1—k)"k dk 
-: “4 


1)) 2 
2 ke2m i 1(] — k2)n 
1 


x {—log k? +-2b(m-+- 1) —2b(m+-4)} dk 


Lh. (v (m+4)) ?0(m-+-1)0 saat 
(Pr (m-+1)} § (m me +-2) 


x {b(m-+-n- bis +-1)—2y(m+4)}. 
= {P(m+4)}P2P(n+1) aad )\* 
(m 


— 1 (m-+-1)I -+-n-+-2) \(n+3) Dye 
mW 

* Quart. J. of Math. 45 (1914), 141-58. 

+ For proofs of this formula, see Watson, Journal London Math. Soc. 4 
(1929), 82-6; Darling, ibid. 5 (1930), 8-9. 


y, 
ot 
7 ( 


m ) 


Since 
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we deduce that 


1 T(n+3) < {T'\(m+4)}* 
G,, =-l 1+. *3 2 
n= OB (n+ P+) 2, Pom-+1)r(m-+n+2) 
x {yo(m-+-n++2)—log(n+ 1)+-yp(m+1)—2p(m+ 4)}. 
Now let M be any fixed integer, not less than 2; it is easy to 
verify that, when m > M, 


0<d(m+n+2)+p(m+1) < 2m+n+3 < 2(m+n+]), 








and 
0 < log(n+1)+ 2(m+4) < (n+1)+2(m+1) < 2(m+n+1). 
Hence 
(P(n+3)}? (P(m+4)}? 
mE (n+), Dr P(m+1)P(m-+n+2) * 
x (lm-+n-+2) —log(n+1) + Hm-+ )—20(m-+4)} 





= Mf 


— AT (n+)? S {M(m+3)} 
mT\(n-+1) D'\(m+1)0(m+n-+1) 


; 2{T'(n-+3)}? — — {T(m+h)P 
71\(n+-1)(M-+-2)(M-+3).. (Mn)  P(m+1)P(m+2) 





m=M 





Pot DPTM+2) wy 4.9.1) —0( 2 ) 

7D(n+1)0(M+n+1) a, $5251) — "lar 
when » is large, since M is fixed. Consequently 

D(n+3 )}? M-1 {T(m+4)¥ o 
mT (n+ 1) Ly P(m-+1)P (m+n +2)’ 

x {p(m+n+2)—log(n+1) br m+ 1) —2y(m+3)}+Of(n+1)-™}, 
and the asymptotic expansion of G, up to the term containing 
(n+1)?-" is obtainable by substituting the asymptotic expansions 


of the functions which occur in this terminating series. 
The asymptotic expansion, as far as I have calculated it, is 





G, , = + log(n+ n+! 


~ 


1 1 ° 
, eile i—~ ———+..1, (8 
G, = |lou(n+1)-+7-+4 log 2 Hn+1)* 19%n+1)* (8) 


and y-+4log 2 = 3-34980 43871. 
4. We now revert to the problem of proving (2) and (3). As a pre- 
liminary we shall show that 


(2n+2)%(pyn43—PynI—(2N+1)%(Pyn—Pyin—») > 0 (W=1,2,8,...). (9) 
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First observe that, when n = 1, the expression on the left in (9) is 
equal to 437—300-+ 964 3 
3a 
and this is positive. It is consequently sufficient to discuss (9) when 
n assumes the values 2, 3, 
It follows from (4) that 


we 
% 


] 
Py = / Ly io—3 fe 


ee {2m(n +1) 1;t) dt 


2m(n+1)+1¥ 


fz —e-i? 


— 


dt, 


2¢ 
so that 


va 
wo . 
16 S ] * sinh mt o 
Pan Pun )j—za — e~m2n 1 dt. 


4m?— 1 sinh t 
1 0 
and hence 


(2n + 2)?( Prenat Pyn)—(22+1)*(pyn—pycn ») 


L 
» 2 -— 
16 * I sinh me «on +4-2)%e m(2n 3 (2n+ 1)2e—mi2n +0) di, 
7 a 4m?—1} sinht | 


mam 0 
It is evident from this result that, in order to prove (9) when 
2, it is adequate to prove that 
a 
* sinh mt ' er 
: {(2n+2)%e-m2n+3"__ (Qn + 1)%e-m2n+0} de > 0 (10) 
sinh¢ 
0 
Now consider (10); for brevity write 
{(2n-+-2)*e m(2n 3” (2n-+-1)? e—m(2n ‘DY cosh(m— 1)é = 
sinh mt 
sinh ¢ cosh(m— 1)t 
For any given values of m and n, it is obvious that ¢(¢) vanishes 
for only one positive value of t, say t= a; and 
d(t)>0, (t<a); d(t)<0, (t>a). 
It is evident that y(t) is a constant when m= 1 or m = 2, and it 
is easy to prove that it is a decreasing function of t when m > 2; for 
_(m—1)sinh 2t—sinh(2m—2)t 


‘(t BBRSE os ws <0, 
x) 2 sinh*t cosh?(m— 1)t 
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when m—1> 1, as may be seen by expanding the numerator in 
powers of ¢t. Hence 

t 
XO 51 <a); (11) 
x(«) 


Now consider [ d(t) dt; we have 
0 
E 2n-+ 2)? 1 
$y ae=' _ epoch Tagaytay A 
one 1 “a 1 
2 |m(2n+1)—(m—1) © m(2n+1)+(m—1) 
_ —m(2n+2)(2n+3) _ m(2n+1)% 
m*(2n+3)?—(m—1)? m?(2n-+-1)?—(m—1)? 
— m{m?(2n+- 1)?(2n+3)—(m—1)*(16n?+ 26n-+ 11)} 
{m?(2n-+-3)?—(m— 1)?}{m?(2n+-1)?—(m—1)*} 
m(m— 1)*{(2n-+ 1)?(2n+-3)—(16n?+ 26n-+ 11)} 
{m?(2n-+ 3)2—(m—1)?}{m2(2n+ 1)?—(m— 1)?} 
i 4m(m— 1)?(n+ 1)(2n?—n—2) 
{m?(2n+3)?—(m—1)?}{m?(2n-+- 1)?—(m—1)*} 
> 0 when n > 2. 














a «o 


Hence, when n > 2, [ d(t) dt > [ {—¢(t)} dt 
0 & 


and a fortiori, by (11), 


t 7 t 
| oo a> [022 
x(2) J 
0 a 
that is to say, [ d(t)x(t) dt > 0. 
0 

Hence (9) is true when n> 2, and it has already been proved 
when n = 1; it is consequently om for all of the specified values of n. 

5. It follows at once from (9) that, when VN > >1, 


Pin—Pn- a, P3(n+1)— Pin 5<: 
1.3...(2n—1))® +i (2n-+ 1)? 
2.4...(2n+2) 
4 


Bo Wicsisies 2n 
P3N—PHN-1)__ _fin—Pwyp = 
<713..2N—1)2 3 < lim 7 QN—1)? a? 
$4... _& WeMERe, 2N 
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so that 4 
(4,,—G,,-1) (n ==}, 2, ae 


7 


Pin Pan-» 


By summing inequalities of the type (12), we get 


4 35 
Pin—Po < —(G,—Go) (n = 1, 2, 3,...), 
7 


and a fortiori 4 
Y € 
Pin <-—G, = Ll, 2, Bch 
: 7 


This completes the proof of (2). 
6. To complete the proof of (3) it is now sufficient to prove that 
a Fan (n = 1, 2, 3....). 
To effect this object, write (9) in the form 
Cus, _. G,.—G,2 
Pi(n4+1)— Pan ~ Pin—PHn-» 
so that 
¢. 
Pint) _] Pin _y 1 — PHm— Pin __] 
Pan P4(n-1) , Phn Pi(n-D 
By summing inequalities of this type, we get 
C.—G, G,—G 


n 


Y Y Y Y Y Y Y 
1—G, _ G, —G,,. ae G,,—G, = G,.—G,_, 7 G BES | 
~ — Un n=i° 





2£<G,—-G, (n = 1, 2, 3....), 





Punt) __} 8 
Pin Po 
that is to say, 
G.. 14—G, —G,, 25—4_¢.. ae 0, 
Pxnti)__ Pr_y 16—47 
Pin Po 
and therefore 
G,44.—G, <@, (PH ae (n = 1, 2, 3,...), 
\ Pin 
which is equivalent to the desired result. 
It may be remarked that the last inequality shows that 
1 rn . 1 4 4 
Pin < Awe See < me Pan =—, 
Gi, Gi, +1 Nw Vn 7 
and this procedure gives an alternative method of proving the result 


of § 5. 
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